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ABSTRACT 


This  stmiy  presents  an  analysis  of  the  fixed  biases  and  variable 
errors  affecting  gun  fire.  The  emphasis  is  placed  on  automatic  cannons 
being  developed  for  use  as  main  guns  for  armored  vehicles. 

It  is  proposed  that  barrel  bend  has  been  overlooked  as  a primary 
cause  of  inaccurate  fire  and  that  automatic  cannons  will  be  particularly 
susceptible  to  this  condition.  A formula  is  derived  for  approximating 
the  magnitude  of  bend  that  can  occur. 

A model  is  presented  for  estimating  miss  distance  in  the  plane  of 
the  target  as  a result  of  fixed  biases  (including  barrel  bend)  and  other 
variable  errors.  A parametric  analysis  demonstrates  the  magnitude  of 
resxilting  miss  distance . Methods  of  avoiding  the  effects  of  barrel 
bend  are  presented. 
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I . INTRODUCTION 
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With  the  introduction  of  vastly  improved  threat  weapons  and  observed 
numerical  imbalances,  the  Department  of  Defense  over  the  last  decade 
initiated  many  programs  to  analyze  existing  and  advanced  state  of  the 
2urt  technologies  for  armored  combat  vehicles  of  the  future . According 
to  the  United  States  Army  Armor  and  Engineer  Board  CQ,  these  programs 
were  designed  to  improve  the  firepower,  mobility/agility  and  protection 
of  future  combat  vehicles  while  concurrently  reducing  weight  and  cost. 

As  noted,  an  integral  feature  of  these  technological  programs  is  the 
development  of  technology  for  improved  firepower  ajid  current  interest  is 
vested  in  the  medium  caliber  automatic  gun.  The  Department  of  Defense 
ainalysts,  especially  research  and  development  analysts,  as  well  as 
tactical  and  strategic  plajiners  must  evaluate  data  on  the  weapon  per- 
formance. This  need  for  methods  of  evaluating  weapons  systems  becomes 
paramoiint  since  our  future  depends  on  intelligent  and  informed  decisions 
at  the  earliest  stages  of  any  program. 

A current  example  began  early  in  the  1970 's  as  the  Defense  Advance 
Research  Projects  Agency  (DARPA)  developed  a program  to  demonstrate  new 
technologies  for  futiare  tanks.  One  of  the  major  component  technology 
areas  in  the  program  is  a vehicle  characterized  by  high  mobility  and 
agility  and  hence  the  program  named,  HIMAG.  A critical  element  in  the 
HIMAG  program  is  the  demonstration  of  the  feasibility  of  a medium  cali- 
ber automatic  gun.  The  weapon  pod  proposed  for  evaluation  is  one 
containing  the  ARES  75  mm,  solid  propellant,  rapid  fire  cannon  using 
APFSDS  and  HE  inert  rounds.  This  75mm  gun  is  designed  to  fire  in 
various  modes  and  inco2T)orates  the  capability  for  both  rapid  or  single 
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shot  fire.  Here  we  have  an  explicit  need  for  a method  to  evaluate  the 
proposed  weapon  system  in  a burst  mode  to  demonstrate  the  technical 
feasibility  of  the  advanced  automatic  cannon. 

However,  Pfleger  and  Cerra;.o  report  that  to  date  there  has  been 
only  limited  effort  directed  towards  providing  reliable,  accucrate  anal- 
ysis of  the  performajice  of  automatic  gun  systems.  Primary  emphasis  has 
traditionally  been  to  evaluate  single  shot  weapon  systems.  Consequently, 
with  the  advent  of  the  automatic  cannon  and  the  ciirrent  need  for  effec- 
tive weapon  analysis  as  generated  by  the  Department  of  Dex'ense,  this 
paper's  objectives  axe  to: 

1 . Develop  a mathematical  model  to  predict  approximate  miss  distance  , 
of  the  projectile  at  the  target  plane  for  medium  caliber  automatic 
cannons  in  order  to  provide  accurate  and  reliable  information  upon 
which  hit  probabilities  can  be  developed  for  analysis  by  DOD. 

2.  Develop  a mathematical  model  to  approximate  barrel  bend  and  use  this 
model  to  determine  the  possible  amount  of  error  caused  by  barrel  bend 
using  the  AEPA  75™  automatic  cannon  as  a test  case. 

3.  Use  the  mathematical  models  developed  in  objectives  "1"  and  "2"  in 
a computerized  parametric  analysis  of  the  ARPA  75™  cannon  under  a 
baseline  M60A1  fire  control  system  approximation  and  the  more  advanced 
XM-1  fire  control  system  approximation. 

A more  detailed  statement  and  presentation  of  the  automatic  cannon 
accuracy  evaluation  problems  are  discussed  in  Chapter  II,  Nature  of  the 
Problem.  The  body  of  the  study  is  presented  in  Chapter  III,  Presenta- 
tion and  Analysis  of  Data.  Chapter  III,  Section  A,  Miss  Distance 
Analysis,  provides  a background  for  the  reader  on  error  sources  in- 
volved with  automatic  cannons  on  combat  vehicles  as  well  as  develops 
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the  mathematical  model  on  miss  distance.  Chapter  III,  Section  B, 

Barrel  Bend  Analysis,  provides  a detailed  analysis  and  mathematical 
model  on  barrel  bend  approximation.  Once  these  mathematical  models 
have  been  formulated,  then  Chapter  IV,  The  Simulation  Model,  employs 
the  models  developed  using  the  ARPA  73im  automatic  cannon  as  a test 
case  under  several  assumed  weather  conditions.  Lastly,  Chapter  V is 
dedicated  to  a summary  and  conclusions  of  the  study  as  generated  during 
I the  research  effort  and  from  analysis  of  the  simulation.  i 


* * : 


II.  NATURE  OF  THE  PROBLEM 


Neaxly  3000  tanks  were  destroyed  during  the  1973  Arab-Israeli  war. 
Post-war  analysis  has  shown  that  the  clear  majority  of  these  tanks  were 
destroyed  by  the  main  gun  of  other  tanks.  Thus,  as  late  as  1973*  the 
tank  remains  the  primary  anti-tank  weapon.  Analysis  also  reveals  the 
total  rounds  fired  considerably  exceeds  the  quantity  expected  to  be 
fired  for  the  number  of  hits  actiially  achieved.  It  can  be  inferred 
from  the  data  that  tank  guns  are  not  nearly  so  accurate  as  commonly 
believed. 

In  the  event  of  a conventional  war,  armored  vehicles  constitute  a 
primary  threat  to  the  secirrity  of  our  allies.  If  the  tank  is  to  be  the 
primary  anti-tank  weapon  in  the  future  or  if  the  tank  is  even  to  re- 
main effective  as  an  anti-tank  weapon,  something  must  be  done  to  im- 
prove the  accuracy  of  tank  gunfire . 

The  primary  effort  toward  improved  accuracy  has  been  the  development 
of  increasingly  complex  and  sophisticated  fire  control  instrumentation. 
Significant  progress  has  been  made  in  minimizing  the  effect  of  some  of 
the  factors  which  contribute  to  miss  distance  at  the  target,  e.g.  range 
estimation . 

A second  approach  to  improving  the  probability  of  hit  is  currently 
being  investigated.  That  is,  to  fire  a burst  of  several  rounds  at  the 
target  instead  of  a single  round.  If  each  round  of  a burst  can  be  fired 
in  precisely  the  same  intended  trajectory,  round-to-round  dispersion 
will  cause  a shotgun  effect  at  the  plane  of  the  target.  The  probability 
of  a hit  in  such  a bxirst  would  clearly  be  higher  than  for  a single 
round . 
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To  demonstrate  this  approach  a weapon  must  be  developed  that  is 
capable  of  burst  fire.  This  means  firing  anti-tank  rounds  at  rates  of 
fire  approaching  120  rounds  per  minute  or,  more  appropriately,  2 rounds 
per  second.  Weapons  of  this  sort  are  being  developed.  The  ARES  75nun 
gun  currently  in  development  has  similar  characteristics.  This  gun  will 
be  used  as  a model  when,  in  the  course  of  this  paper,  it  is  necessary 
to  utilize  the  specific  characteristics  of  a rapid  fire  gun. 

Given  a gun  with  rapid  or  burst  fire  capability,  it  must  still  be 
demonstrated  that  each  round  of  a burst  is  fired  with  the  same  degree  of 
accuracy.  This  requires  the  barrel  to  be  pointing  in  the  same  direction. 
It  must  be  demonstrated  that  the  recoil  or  trunnion  forces  can  be  trans- 
ferred through  the  entire  vehicle  and  the  gun  platform  and  gun  can  be 
returned  to  the  proper  position  in  the  time  available.  Rapid  fire  also 
means  higher  bairrel  temperatures.  High  temperatures  can  affect  barrel 
life  and  accuracy.  Safety  is  another  consideration  in  view  of  the 
possibility  of  ignition  of  the  propellant  due  to  chamber  temperature  or 
"cook  offs."  It  is  clear  that  there  is  much  to  be  considered  and 
demonstrated  in  the  development  of  a rapid  fire  tank  gun. 

A third  approach  is  implied  in  the  foregoing  but  should  be  addressed 
separately.  Determine  the  factors  which  contribute  to  miss  distance 
that  have  not  been  eliminated  or  minimized  to  date  and  find  ways  to 
compensate  for  or  avoid  their  effects.  It  is  felt  that  far  too  little 
attention  is  being  given  to  two  factors  which  make  ’significant  con- 
tribution to  miss  distance,  i.e.  round- to-round  dispersion  and  barrel 
bend. 

Round- to-round  dispersion  is  primarily  the  result  of  differences  in 
the  physical  characteristics  of  each  round  which  ad’fect  the  trajectory 
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of  the  round.  This  dispersion  is  well  known  to  those  familiar  with 
gunfire  in  general.  Unfortunately,  this  error  has  been  deemed  "random" 
which  implies  "unavoidable."  It  is  felt  that  further  efforts  to  mini- 
mize this  factor  are  justified. 

Barrel  bend  is  a little  known  factor.  Few  recognize  that  gun 
baorrels  do,  in  fact,  bend.  Even  fewer  have  considered  the  possible 
magnitvide  of  the  resulting  miss  distance.  It  is  understood  that  weapons 
are  being  developed  by  West  Germany  and  Great  Britain  that  may  avoid 
barrel  bend  or  its  effect.  There  is  no  information  available  to 
indicate  that  such  a weapon  is  being  considered  in  this  country. 

A primary  purpose  of  this  paper  is  to  demonstrate  that  barrel  bend 
can  cause  a significant  error,  particularly  in  view  of  the  higher 
barrel  temperatures  that  will  occur  with  rapid  fire  systems . 
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III.  PRESENTATION  AND  ANALYSIS  OF  DATA 


A.  MISS  DISTANCE  ANALYSIS 
1 . Introduction 

This  section  describes  an  analytical  method  formulated  for  use 
in  the  evaluation  of  the  performance  of  automatic  gun  systems  for  combat 
vehicles  against  stationary  ground  taurgets.  The  need  for  such  analytical 
effort  is  readily  apparent,  especially  in  view  of  the  current  emphasis 
on  armored  combat  vehicle  technology  and  the  HIMAG  program  evaluating 
the  use  of  the  75n™  automatic  gun. 

The  first  step  in  the  analysis  was  to  determine  a significant 
effectiveness  measure  for  the  automatic  gun.  Traditionally,  hit  proba- 
bility (Pj^)  has  been  used  in  this  case  and  considerable  effort  in 
measuring  this  effectiveness  has  been  accomplished  and  guidelines 
established  for  first  round  or  single  shot  hit  probability  by  Frankford 
Arsenal  D.-a-  However,  with  the  advent  of  the  automatic  gun  which  can 
fire  in  a burst  or  rapid  fire  mode,  the  hit  probability  of  subsequent 
rounds  fired  in  rapid  succession  must  be  addressed. 

Secondly,  it  is  felt  that  any  analysis  in  this  area  must  begin 
by  determining  the  impact  point  of  each  round  fired,  either  by  single 
shot  or  burst  fire.  By  mathematically  modeling  the  impact  point  of  each 
round  fired  at  the  target  plane  in  relation  to  the  center  of  the  target, 
a base  would  be  provided  for  any  additional  analysis  in  hit  probability 
evaluations.  Logically,  the  farther  the  impact  point  of  a round  from 
the  center  of  target,  the  lower  the  overall  hit  probability.  The  use 
of  all  error  sources  is  designed  to  simulate  the  conditions  which  can 
be  expected  under  a variety  of  actual  combat  conditions.  According  to 
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Brcxikln  ^33*  this  measure,  Quasi-Combat  Condition,  is  Intended  to  reflect 
as  realistically  ais  possible  the  combat  potential  of  the  system  in  re- 
gard to  delivery  accuracy. 

Lastly,  under  the  scenario  input  to  the  simulation  in  Chapter  IV, 
two  terms  mentioned  by  Cerrato  D]  necessitate  reiteration.  An  "occasion" 
is  one  tactical  engagement  between  vehicles  or  stationary  implacement. 

A "burst"  is  a series  of  rounds  fired  in  rapid  succession  with  the  same 
point  of  aim. 

2.  Classification  of  Error  Types 

Since  the  standard  target  is  two-dimensional,  the  effects  of 
both  horizontal  (X)  and  vertical  (y)  error  sources  must  be  considered. 
These  gunnery  error  sources  cause  certain  errors  in  different  magnitudes 
and  in  both  directions  at  the  target  plane.  Pfleger  and  Bibbero  re- 
port that  this  leads  to  the  assumption  that  the  horizontal  and  vertical 
errors  are  independent  of  each  other.  This  method  is  valid  if  the 
target  is  aligned  with  the  coordinates  of  the  plane  in  which  the  errors 
are  considered. 

From  the  classification  of  error  types  presented  in  Frankford 
Arsenal  ^,5^  » gunnery  errors  can  be  classified  into  three  types  with 
respect  to  automatic  cannons  on  combat  vehicles. 

a.  Fixed  Biases  are  a consequence  of  system  design  and  use,  and 
are  constant  for  a given  range.  In  the  fire  control  systems  in  use 
today,  correction  is  made  for  many  of  these  fixed  biases. 

b.  Variable  Biases  are  considered  to  be  constant  for  a specific 
engagement  and/or  burst  of  fire.  They  vary  in  a random  fashion  but  not 
from  round  to  round.  Included  in  this  category  ares 

(1)  Occasion- to-Occasion  biases  are  errors  which  change  over 
time,  but  change  is  so  slow  that  they  can  be  considered  a constant  value 
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over  the  period  of  an  engagement.  Generally,  these  errors  are  induced 
by  the  weapon  system  design  a«id  the  tactical  situation. 

(2)  Burst- to-Burst  biases  are  considered  to  have  different 
values  for  e?.ch  burst  fired  during  the  engagement  and  are  considered  to 
vary  in  a random  manner. 

c.  Random  Errors  vary  from  round  to  round  in  such  an  unpredic- 
table manner  that  no  correction  can  be  introduced  by  the  fire  control 
system.  These  random  errors  are  classified  as  Round- to-Round  errors  and 
are  the  result  of  variations  in  ammunition  rounds  and  characteristics  of 
the  weapon  system. 

As  can  be  seen  by  the  error  classification,  the  mean  (p-)  of  the  dis- 
tribution of  any  round  fired  is  determined  by  the  fixed  biases  of  the 
system.  That  is,  the  system  bias  offsets  or  shifts  the  position  of  the 
stauidard  normal  curve  from  the  target  center  of  mass  to  the  center  of 
impact,  (ps)  The  variable  biases  which  are  considered  constant  for  the 
particular  engagement  and  each  burst  within  the  engagement  cause  a dis- 
persion (O^)  about  the  impact  point.  The  random  errors  provide  an  index 
of  the  additional  dispersion  (o^)  distance  from  the  meain  of  the  dis- 
tribution of  any  round  fired  in  the  biirst. 

In  the  firing  of  an  automatic  cannon,  it  is  compulsory  to  calculate 
the  effects  of  several  of  the  variable  biases  which  have  traditionally 
been  listed  only  as  a system  variable  bias.  BrodJcLn  02  reported  that 
if  the  calculated  effects  of  these  variable  biases  can  be  accurately 
measiired  that  this  provides  a more  representative  impact  point.  This 
is  extremely  important  in  automatic  cannon  gunnery,  since  deflection  of 
the  muzzle  and  muzzle  velocity  changes  caui  realistically  be  encountered 


in  the  burst  of  fire.  As  a result,  if  several  variable  biases  are  cal- 
culated they  must  be  Included  as  an  additional  fixed  bias  of  the  system 
and  deleted  as  a source  of  variable  bias . 

3.  Miss  Distance 

The  point  of  impact  at  the  target  plane  can  be  addressed  by 
using  the  X and  Y axes  for  reference.  In  most  cases  the  gunner  vdll 
have  a line  of  sight  to  the  target  at  the  center  of  mass  of  the  target. 
This  point  will  be  at  coordinates  0,0  of  the  target  plane.  The  offset 
from  the  target  centerpoint  when  the  projectile  arrives  at  the  target 
plane  is  due  to  fixed  system  biases  and  will  be  defined  as  the  expected 
”Mlss  Distance."  This  expected  miss  distance  is  defined  in  both  the 
azimuth  (p-X)  and  elevation  (p*!)  directions  at  the  target.  Synony- 
mously this  miss  distance  is  termed  "X-error"  and  "Y-error"  for  azimuth 
and  elevation  respectively.  (See  Figure  1) 


Line  of 
Sight 


Figure  Ij  Miss  Distance  Representation 
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This  study  attempts  to  mathematically  model  the  miss  distance 


for  any  round  at  the  target  plane  by  considering  the  system  fixed  biases, 
as  well  as  accounting  for  any  variable  biases  which  can  be  accurately 
measured  during  the  subject  firing  engagement.  In  order  to  compute 
X-error,  the  total  displacement  of  the  round  from  a given  muzzle  posi- 
tion must  be  determined.  Dickinson  ^63  reports  this  lateral  displace- 
ment from  the  line  of  fire  is  caused  by  crosswind,  rotation  of  the  earth 
and  drift  as  addressed  below. 

a.  Drift  (Xj^)  is  caused  by  aerodynamic  forces  which  act  on  the 
projectile  along  its  intended  route  of  flight.  Most  rounds  which  exhib- 
it drift  are  spin-stabilized  and  have  high  spin  rates.  Any  drift  de- 
flection is  positive  and  arises  from  the  clockwise  spin  of  the  projectile. 

b.  The  rotation  of  the  earth  (Xg.Yg)  caiises  a system  fixed  bias 
as  a result  of  Coriolis  acceleration.  However,  the  error  is  assumed  for 
this  study  to  fall  into  the  category  of  Occasion-to-Occasion  variable 
bias  due  to  random  characteristics.  (C7Xg,0’Yg) 

c.  Any  crosswind  (X^)  force  present  during  the  projectile's 
flight  will  cause  a deviation  in  azimuth  from  its  intended  trajectory. 
Normally  this  displacement  is  listed  as  a variable  bias  but  it  is  felt 
that  this  must  be  addressed  in  the  X-error  computation  for  accuracy  of 
trajectory  under  various  Quasi-Combat  'rfeather  Conditions,  since  weather 
will  be  a definite  factor  influencing  automatic  gun  firing. 

d.  Additionally,  a deflection  of  the  muzzle  (X^j)  will  cause  a 
displacement  of  the  round  in  the  same  direction  at  the  target.  This  de- 
flection is  caused  primarily  by  tube  or  barrel  bend,  improper  gun  align- 
ment, and  trunnion  force  effect  of  previously  fired  rounds  in  the 
burst. 
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In  past  analysis  of  hit  probability,  Jump  (X^)  and  Parallax  (X_) 

jk 

have  been  liswa  as  a fixed  bias  of  the  gun  system.  However  when  moder- 
ately sophisticated  fire  control  systems  are  employed,  these  errors  are 

usually  fully  corrected  according  to  Pfleger  a«id  Bibbero  [j+J.  As  such, 

2 

these  fixed  biases  are  considered  a round  dispersion  (O'  ) factor  and 
will  not  be  a factor  in  any  X-error  model  for  evaluation  of  advanced 
technology  gun  systems. 

The  derivation  of  miss  distance  in  the  Y direction  (Y-error)  from 
the  line  of  fire  is  caused  by  the  following  factors  reported  by  Pfleger 
and  Bibbero  [fO  • 

a.  Deviation  from  projectile  standard  muzzle  velocity  (Y^)  is 
the  result  of  powder  temperature  variations,  tube  heat  transfer  and 


tube  bend. 

b.  Any  deviation  from  the  standard  air  density  (Y^)  will  cause 
drag  effects  of  the  air  mass  through  which  the  projectile  passes  toward 
the  target.  Accordingly,  any  percent  change  from  standard  air  density 
due  to  meteorological  effects  will*^ncrease  or  decrease  the  range  of  the 
projectile . 

c.  Deviation  from  standard  air  temperature  (Y^)  will  cause  an 
increase  or  decrease  in  range  of  projectile  similar  to  air  density 


changes. 

d.  Any  rangewind  (Y^)  fqgce  will  cause  a deviation  in  elevation 
from  the  projectile's  nominal  trajectory. 

e.  Range  error  fixed  bias  (Y^)  is  caused  by  inaccurate  range 
estimation.  Pfleger  and  Bibbero  report  that  for  the  laser  range- 
finder this  may  be  taken  as  zero,  but  when  distance  measurements  are 
made  by  the  human  eye,  the  fixed  bias  cannot  be  stated  with  certainty. 
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Thus,  range  estimation  should  only  be  considered  as  a variable  bias 


(CTYp)  since  the  crew's  estimation  of  the  range  will  probably  change  in 
a random  fashion  for  each  engagement. 

f.  Any  elevation  of  the  muzzle  (Y^^)  deviating  from  the  appro- 
priate superelevation  setting  for  the  specific  range  will  increase  or 
decrease  the  projectile  rainge.  This  deviation  is  caused  by  gunner 

I 

superelevation  error,  barrel  bend,  and  trunnion  force  effect  from  previ-  ; 

ously  fired  rounds  in  the  burst  of  fire.  ! 

As  was  previously  noted  in  the  X-error  biases,  there  are  certain 
biases  which  have  traditionally  been  considered  fixed  but  due  to  modem 
fire  control  systems,  the  assumption  is  made  that  they  are  fully  coirec- 
ted.  Likewise,  in  computing  Y-error  at  the  target,  the  fixed  biases  of 
Jump  (Yj)  and  Parallax  (Yp)  are  considered  to  be  zero  and  thus  fully 
corrected. 

In  this  analysis  the  use  of  X-error  and  Y-error,  miss  distances, 
give  the  amalyst  a good  indication  of  the  point  of  impact  for  either  a 
single  shot  or  a burst  of  several  rounds.  Chapter  m.  Section  B in- 
cludes a discussion  of  the  impact  of  barrel  bend  on  the  automatic  gun 
and  the  resultant  position  of  the  muzzle  at  time  of  firing,  at  which 
time  the  approximate  point  of  impact  cam  be  determined  for  each  round. 

By  knowing  the  mean  ((-t)  of  the  standard  normal  curve  determined  by  the 
fixed  bias  of  the  system  and  the  dispersion  about  that  point  in  the  form 
of  a variance  (<J  ) , then  integration  of  an  assumed  distribution  over 
the  target  area  is  required  to  generate  hit  probabilities. 

4.  Mathematical  Model  for  Miss  Distance 

The  computation  of  miss  distance  requires  a mathematical  model 
which  will  express  miss  distance  as  a function  of  the  sources  of  error 
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^ 

affecting  the  firing  trajectory  of  the  projectile . The  variables  which 
have  an  effect  on  miss  distance  at  the  target  plane  have  previously 
been  enumerated  in  this  chapter.  Using  these  variables  in  conjunction 
with  the  appropriate  firing  tables  for  a particular  projectile,  the  model 
can  be  derived  to  give  a good  approximation  of  miss  distance  at  the 
particular  range  in  q,uestion. 

This  model  is  based  on  the  following  assumptions: 

a.  Variable  biases  can  be  combined  with  the  fixed  system  biases 
to  refine  impact  point  (p)  if  the  value  can  be  measured  directly  and 
hence  the  effect  determined. 

b.  Jump  (X^,Y^)  fixed  bias  equals  zero  due  to  modem  fire  con- 

J J 

2 

trol  systems  and  is  treated  only  as  an  element  of  system  dispersion  (CT  ) . 

c.  Paxallax  (Xp,Yp)  fixed  bias  equals  zero  due  to  modem  fire 
control  systems. 

d.  Range  Estimation  (Y^)  fixed  bias  equals  zero  and  is  addressed 
only  as  an  element  of  the  system  dispersion 

e . The  firing  table  data  used  is  accurate  for  the  particular 
round  in  question  and  will  provide  appropriate  variable  coefficients. 

f.  The  horizontal  (X)  and  the  vertical  (Y)  errors  are  indepen- 
dent of  each  other. 

g.  Fire  control  residual  sensing  error  is  accurately  assessed. 

h.  Deflection  of  the  muzzle  as  a result  of  trunnion  forces 
equals  zero  for  any  round  fired.  Tills  is  in  fact  not  realistic  but 
little  research  in  this  area  is  available  for  the  study. 

The  miss  distance  model  for  X-error  (X^)  is  developed  by  con- 
sidering the  variables  which  affect  the  lateral  displacement  of  the 
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round  from  the  line  of  fire.  The  vaxiables  previously  listed  in  Chapter 


III  au^ei 


a.  Drift 


b .  Crosswind 


c.  Muzzle  lateral  displacement  resulting  from  inaccurate  gunner 
laying,  barrel  bend,  and  effect  of  trunnion  forces. 

The  predicted  point  of  impact  of  the  projectile  in  meters  at  the 
target  can  be  modeled  using  these  variables  in  conjunction  with  firing 
table  data  for  the  round  and  weapon  system  as  illustrated  below. 


^T  ' °D  ^2^ 

where  the  coefficients  a.  and  a^  are  given  by: 

X At 

^ d Drift 

1 dm/s  crosswind 


which  is  the  change  in  projectile  drift  for  a unit  change  in  crosswind 
in  meters  per  second  from  the  firing  table 

^ _ d Lateral  displacement 
2 ” d Lateral  muzzle  position 

which  is  the  change  in  lateral  displacement  of  projectile  for  a unit 
change  in  lateral  muzzle  position  expressed  in  mils. 

The  coefficient,  a^,  is  computed  by  use  of  a simple  geometric 
triangulation  for  direct  fire  weapons  as  shown  in  Figure  2.  The  width 
of  the  change  can  be  computed  by  the  following  formula. 


Width  (a.)  = MILS  X 


Range 


For  example,  at  1000  meters,  a 1 mil  (^)  change  in  muzzle  angle  toward 
the  target  will  result  in  a 1 meter  lateral  displacement  of  the  round 
at  the  target  plane . 
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Figure  Zt  Geometry  of  Muzzle  Lateral  Displacement 

The  constant  and  variable  derivations  are: 

Gjj  = Drift  of  the  round  from  firing  table  data  at  the  desired 

range 

X,^  = Grosswind  in  meters  per  second 

(1)  if  crosswind  is  left  to  right  = positive 

(2)  if  crosswind  is  right  to  left  = negative 

= Muzzle  lateral  displacement  in  mils  resulting  from  the 
summation  of  gunner  error  (X^) , barrel  bend  (Xg) , and  trunnion  force 
deflection  (Xp) . Thus: 

* ’S  * ’'p  (3) 

This  model  for  X-error  is  valid  only  if  no  fire  control  system 
wind  sensors  are  employed.  If  sensors  are  employed  then  the  model  must 
be  modified  to  account  for  the  accuracy  of  the  sensor  since  any  correc- 
tion is  a function  of  sensor  and  correction  accuracy.  The  X-error  will 
be  the  residual  error  remaining  after  sensing  and  corrective  action  to 
the  system. 
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Past  data  on  wind  magnitudes  indicates  that  an  average  of  approx- 
imately 50^  wind  sensing  and  correction  is  possible  but  this  is  again  a 
function  of  the  sensor  accuracy.  Thus  for  the  system  being  evaluated, 
the  sensing  error  miist  be  calculated  correctly  and  applied  to  Equation  2 
in  the  following  manner. 


where  S„  is  the  approximate  residual  sensing  error  for  crosswind. 

n 

The  miss  distance  model  for  Y-error  (Y^)  is  developed  by  con- 
sidering the  variables  which  affect  the  vertical  displacement  of  the 
round  from  the  line  of  fire  to  center  of  target.  These  variables  were 
previously  listed  in  Chapter  III  as: 

a.  Deviation  from  projectile  standard  muzzle  velocity 

b.  Deviation  from  standard  air  density 

c.  Deviation  from  standard  air  temperature 

d.  Range wind 

e.  Muzzle  vertical  displacement  resulting  from  inaccurate 
gunner  laying  and  barrel  bend. 

The  model  below  can  be  used  to  estimate  the  Y-error  in  meters  at  the 
target  plane  by  inserting  these  variables  into  the  general  equation: 


Yt  - -C^  + a^Yy  + 


^2^D  ^ ^3'^A  * 
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which  is  the  change  in  elevation  for  a unit  change  in  meters  per  second 
from  standard  muzzle  velocity  as  published  for  the  projectile. 

d elevation 

^2  “ d percent  standard  air  density 

which  is  the  change  in  elevation  for  a 1%  change  from  standard  air 
density  as  listed  in  the  firing  tables 

d elevation 

a s — ■ I 

3 d percent  standard  air  temp 

which  is  the  change  in  elevation  for  a 1%  change  from  standard  air 
temperature  as  listed  in  the  firing  tables 

d elevation 

“ d meter/ sec.  rangewind 

which  is  the  change  in  elevation  for  a unit  change  of  rangewind  as 
listed  in  the  firing  tables 

d elevation 
5 d muzzle  elevation 

which  is  the  change  in  elevation  for  a 1 mil  elevation  of  the  muzzle 
per  firing  tables . 

Once  the  general  model  is  established  for  the  specified  round  and  range, 
then  Equation  5 is  used  to  establish  Y-error  at  the  target  plane  by  in- 
put to  the  following  variables; 

Yy  = the  change  in  meters  per  second  from  standard  muzzle 
velocity  of  the  projectile  specification  data 

(1)  if  the  muzzle  velocity  increases,  then  Y^  is  positive 

(2)  if  the  muzzle  velocity  decreases,  then  Y^  is  negative 
Yjj  = the  percentage  increase/de crease  from  the  standard  air 

density  used  in  computing  firing  table  data  since  firing  tables  are 


( 

1 

i 
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computed  using  the  International  Civil  Aviation  Organization  (iCAO) 
standard  atmosphere . Standatrd  density  must  be  adjusted  accordingly  as 
altitude  above  sea  level  and  geographic  area  are  changed.  The  model  for 
calculating  the  air  density  in  gms/cc  is  defined  in  Humphreys  QtJ  as: 

_ 1.2930  * 10'^ * B-0.378e 

1 + 2.039  * 10"^(t-32)  29.921  (6) 

where 

t = air  temperature  in  degrees  Farenheit 
B = total  pressure  in  inches  of  mercury 
e = vapor  pressure  in  mm  of  mercury 

By  use  of  the  tables  in  Smithsonian  , a vapor  pressure  can  be  deter- 
mined for  any  temperature  and  relative  humidity  (RH) . Similarly,  B 
(the  total  pressure  in  inches  of  mercury)  can  be  determined  at  any 
standard  atmospheric  altitude  from  Humphreys  . 

= the  percentage  increase/decrease  from  the  standard  air 
temperature.  The  ICAO  standard  air  temperature  at  sea  level  (O  altitude) 
is  15°G  and  must  be  adjusted  for  altitude  variation. 

= rangewind  in  meters  per  second 

(1)  if  the  direction  of  the  wind  is  from  the  target  toward 
the  combat  vehicle  along  the  line  of  fire  then  Y.^  is  negative 

(2)  if  the  direction  of  the  wind  is  from  the  combat  vehicle 
toward  the  target  along  the  line  of  fire  then  Y,^  is  positive . 

Yj^  = the  verticle  muzzle  displacement  in  mils  (;li)  resulting  from 

gunner  superelevation  adjustment  (Y_),  barrel  bend  (Y_) , and  trunnion 

force  elevation  (Y  ) . Thus, 

F 


As  with  the  mathematical  model  for  computing  miss  distance  in 
the  X-direction,  this  model  similarly  is  valid  only  when  no  sensors  are 
used  for  sensing  wind,  air  temperature  and  air  density.  When  sensors 
are  used,  the  residual  error  is  a function  of  the  sensing  accuracy  of 
the  fire  control  system.  For  Incorporation  of  sensors  into  the  model, 
the  percentage  of  each  residual  error  source  must  be  calculated.  After 
calculation  of  these  percentage  values  for  the  variables  in  question, 
the  Y-error  model  becomest 

Yt  = -C^  + a^Y^  + a^CYjj  * S^)  + a^(Y^  * S^)  + 

where  S^,  and  aire  the  percentage  of  residual  sensing  error  for 
air  density,  air  temperature,  and  rangewind  respectively. 

5.  Rovuid  Dispersion 

By  using  the  mathematical  models  previously  developed  for  both 
X and  Y Miss  Distance,  the  theoretic  impact  point  at  the  target  plane 
can  be  determined  for  any  round  fired  either  in  a single  shot  or  burst 
fire  mode.  The  X and  Y-error  constitute  an  offset  between  the  mean  of 
the  normal  curve  (center  of  impact)  and  the  center  of  the  target 
according  to  Pfleger  and  Bibbero  H . It  can  be  seen  that  the  greater 
the  expected  miss  distance  due  to  fixed  biases,  the  lower  the  hit 
probability. 

In  Chapter  IH,  Section  A2  all  classes  of  gunner/fire  control 
system  errors  were  defined.  The  magnitude  of  these  errors  - both  fixed 
bias  or  miss  distance,  and  variable  errors  - is  affected  by  the  sophis- 
tication of  the  fire  control  system  used. 
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Since  sophistication  of  the  basic  miss  distance  models  has 
already  been  accomplished  to  account  for  fire  control  system  product  im- 
provement, the  next  area  of  concern  is  the  total  variation  of  the  system 
about  the  target  centerpoint  of  each  round  fired.  Here  aigain,  it  is  felt 
that  fire  control  system  product  improvement  must  be  addressed. 

The  ARPA  HIMAG  program  is  working  with  several  baseline  fire 
control  system  configurations  to  study  the  effects  of  component  improve- 
ments cind  in  defining  minimum  performance  requirements.  For  purposes  of 
this  analysis  it  was  felt  that  the  model  could  be  exercised  by  using  the 
following  two  HIMAG  baseline  configurations  as  listed  in  General  Motors 
Corporation  [al 

a.  M60A1  Approximationj 

Slaved  Ballistic  Reticle  Az,El 
Superelevation  Only  - No  Sensors 
Manual  Range  Estimations 
Direct  View  Day  Only  OPTICS 

b.  XM-1  Approximation: 

Independently  Stabilized  Gun  Director  Az,El 

Full  Solution  Gun  Offsets 

Laser  Rangefinder 

Sensor  and  Manual  Estimates 

Direct  View  Optics  (Day) , FIR  Biocular 

Using  the  classes  of  error  from  Chapter  III,  Section  2 auid  the 
theory  of  error  components  considered  i n Pfleger  and  Bibbero  , a 
breakdown  of  these  errors  by  class  can  be  made . 

"Round- to-Round"  errors  (which  axe  totally  random  since  no 
correction  can  be  made  to  the  system  by  the  fire  control  system)  are 


caused  byt 

a.  Ammo  Dispersion  (crX^iCY^) 

b.  Jump  Dispersion  (<JXj,(JYj) 

These  rajidom  errors  will  be  the  same  for  any  type  fire  control  system 
and  range,  Round-to-round  errors  take  on  different  values  for  each 
round  fired  in  a burst  and  Cerrato  [5]  considers  them  to  be  independent 
and  normally  distributed.  The  combined  standard  deviation  in  one  dimen- 
sion for  these  random  errors  is  calculated  from  the  component  error 
sources  where 

and 


"Burst- to-Burst"  biases  aire  random  and  by  their  very  nature  will 
take  on  a constant  value  for  the  entire  burst  of  fire.  These  errors  are 
caused  by: 

a.  Gun  Pointing  (Laying)  (CT  X^^,  CY^^) 

b.  Visual  Resolution  (crX^tCTY^) 

It  is  estimated  by  Cerrato  bl  that  these  errors  will  be  independent  ajid 
normally  distributed  and  the  combined  standard  deviations  are  of  the 
form: 

+ (11) 

and 


"Occasion-to-Occasion"  biases  chajige  so  slowly  that  they  can  be 
considered  constant  for  the  period  of  engagement  against  a stationary 
target.  Falling  into  this  category  arej 

a.  Boresight  Loss  (crX^.CTYg) 

b.  Fire  Control  (O'XpjCrYp) 

c.  Ranging  (CTY^^) 

d.  Cant  (<XXj,,(TY^) 

e.  Angle  of  Site  (cT'Yg) 

f.  Earth  Rate  (crX^jCTYg) 

g.  Tube  Droop  (cTY^j) 

h.  Zeroing  (1200  M)  (aX2,crY2) 

For  approximation  purposes  these  errors  can  be  considered  independent 
amd  normally  distributed,  and  the  combined  standard  deviations  are  of 


the  formt 


<^*00  = * 


These  errors  are  not  presented  as  an  all  inclusive  list  for 
every  weapon  system  in  use.  They  represent  the  authors'  best  list  for 
use  in  approximating  the  fire  control  systems  being  presented  in  this 
study.  Additionally,  the  computations  performed  in  computing  the  in- 
dividual standard  deviation  values  are  based  on  guidelines  established 
in  Pfleger  and  Bibbero  M and  are  approximations  only.  Accordingly, 
more  exact  values  should  be  used  and  additional  errors  accounted  for 
when  available  for  use. 
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After  categorizing  the  various  errors  for  analysis,  an  approxi- 
mation of  the  system  error  magnitudes  for  total  dispersion  about  the  miss 
distance  point  can  be  computed.  Tables  I and  II  were  generated  by 
combining  assumed  firing  table  data  with  suggested  error  magnitude 
computations  in  Pfleger  and  Bibbero  W for  both  the  ARPA  APFSDS  and 
HE  rounds. 

6 . Model  Application 

Using  assumed  firing  table  data  for  both  the  ARPA  APFSDS  and  HE 
rounds  listed  in  Appendix  A,  the  miss  distance  model  can  be  applied  to 
both  the  m60A1  and  XM-1  fire  control  system  configurations.  By  gener- 
ating the  constant  and  coefficients  for  the  model,  these  equations  can 
then  be  used  in  the  simulation  of  miss  distance  presented  in  Chapter 
IV. 

To  compute  X-error,  the  mathematical  models  developed  in 
Chapter  III,  Section  A4  are  used.  (See  Equation  1 for  M60  FCS  con- 
figuration and  Equation  4 for  XM-1  FCS  approximation) . Inserting 
appropriate  firing  table  data  at  the  desired  range  into  the  model  re- 
sults in  Miss  Distance  formulas  for  use  in  simulation.  These  formulas 
are  summarized  in  Table  III  and  IV. 

The  mathematical  models  previously  developed  in  Chapter  UI, 
Section  A4,  are  used  for  computing  Y-error  at  the  target  plane.  (See 
Equation  5 for  M60  FCS  configuration  and  Equation  8 for  XM-1  FCS  approxi- 
mation) . Tables  V aind  VI  were  generated  for  the  ARPA  APFSDS  and  HE  pro- 
jectiles by  inserting  assumed  firing  table  data  into  the  basic  mathe- 
matical models.  Additionally  it  is  assumed  that  residual  sensing  errors 
are  1%,  2%  and  3^  respectively  for  density,  air  temperature,  and 
range wind . 
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TABLE  I 

FIRE  CONTROL  SYSTEM  ERROR  MAGNITUDES 
APFSDS  ROUND 
QUASI-COMBAT  CONDITIONS 


M60A1  APPROXIMATION 

XM-1  APPROXIMATION 

1500  M 

2500  M 

1500  M 

2500  M 

Error  Source 

Azimutl 

Elev 

Azimutl 

Elev 

Azimuth 

Elev 

Azimuth 

Elev  ■ 

Round-to- 

Round 

.320 

.280 

.320 

.280 

.320 

.280 

.320 

.280 

.220 

.300 

.220 

.300 

.220 



.300 

.220 

.300 

.254 

.254 

•153 

.153 

.020 

.070 

.020 

.070 

.060 

.060 

.060 

.060 

.060 

.060 

.060 

.060 

Dccasion-to- 

Occasion 

.080 

.080 

.080 

.080 

.080 

.080 

.080 

.080 

.070 

.070 

.070 

.070 

.070 

.070 

.070 

.070 

- 

.624 

- 

1.120 

- 

.006 

- 

.006 

.252 

.001 

.436 

.001 

- 

.011 

.019 

.011 

.075 

.075 

.125 

.125 

.075 

.075 

.125 

.125 

.212 

.220 

.294 

.242 

.170 

.185 

.198 

.210 

NOTES I 


(1)  All  error  values  represent  one  standard  deviation  (1  ) of  the  error 

source  measured  in  mils 


(2) 


Y X „ 1000 


3 


1000 


where  a^  and  ^ target  size  in  each  dimension  (mils) 

R = Target  Range  (meters) 

(3)  Gun  laying  errors  in  XM-1  FCS  configuration  reflect  testing  results 
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TABIE  II 

FIRE  CONTROL  SYSTEM  ERROR  MAGNITUDES 
HE  ROUND 

QUASI-COMBAT  CONDIHONS 


m6oai  approximation 


XM-1  APPROXIMATION 


1500  M 


2500  M 


1500  M 


2500  M 


Elev  lAziiiiutn  Elev  lAzirauthl  Elev  I Azimuth!  Elev 


.320 

.280 

.320 

.280 

.320 

.280 

.320 

.280 

.220 

.300 

.220 

.300 

.220 

.300 

.220 

.300 

.254 

.254 

.153 

.153 

.020 

.070 

.020 

.070 

.060 

.060 

.060 

.060 

.060 

.060 

.060 

.060 

.080 

.080 

.080 

.070 

.070 

.070 

1.684 

- 

3.792 

.001 

1.144 

.001 

.025 

- 

.049 

.075 

.125 

.125 

.220 

.605 

.242 

NOTES: 

(1)  All  error  values  represent  one  standard  deviation  (1  ) of  the  error 

source  measured  in  mils 

U)  x,=  ^.  1220  ^,  = ^*^252 

where  aj^  and  a^  = -I  tairget  size  in  each  dimension  (mils) 

R = Target  Range  (meters) 

(3)  Gun  laying  errors  in  XM-1  FCS  configuration  reflect  testing  results 
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TABLE  III 

MISS  DISTANCE  X-ERROR  FORMULAS 
APFSDS  ROUND  UNDER  QUASI -COMBAT  CONDITIONS 


TABLE  IV 

MISS  DISTANCE  X-ERROR  FORMULA 
HE  ROUND  UNDER  QUASI-COMBAT  CONDITIONS 

M60A1  FCS  APPROXIMATION 
1500  Meters 

X^  = 0 + .197Xy  + 1.5Xj, 

2500  Meters 

X^  = 0 + .6li^Xy  + 2.5\ 

XM-1  FCS  APPROXIMATION 
1500  Meters 

x^  = 0 + .197(Xy  * .03)  + 1.5X^ 

2500  Meters 

X^  = 0 + .614(X^  ♦ .03)  + 2.5\ 

NOTEj  Residual  Sensing  Error  for  Crosswind  Assumed  = Jfo 
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TABLE  V 

mss  DISTANCE  Y-ERROR  FORMULA 
APFSDS  ROUND  UNDER  QUASI-COMBAT  CONDITIONS 

m6oai  fcs  approximation 
1500  Meters 

Y^=  -4.179  + .006Y^+  (-.0025Yjj)  + (-.0015Y^)  + 

.OOOY^^  + 1.473Y„ 

2500  Meters 

Y^  = -12.007  + .0144Y^  + (-.OllY^^)  + (-.003Y^)  + 
.0015Yy  + 2.454Yj^ 

XM-1  FCS  APPROXIMATION 
1500  Meters 

Y^  = -4.179  + .006Y^  + -.0025(Yjj  + .01)  + 

-.0015(Y^+  .02)  + .000(Y,^  ♦ .03)  + 1.^73Yj^ 

2500  Meters 

Y^  = -12.007  + .0l44Yy  + -.01l(Yjj  * .01)  + 

-.003(Y^  + .02)  + .0015(Y^  ♦ .03)  + 2.454Yj^ 


TABLE  VI 

MISS  DISTANCE  Y -ERROR  FORMULA 
HE  ROUND  UNDER  QUASI -COMBAT  CONDITIONS 


M60A1  FGS  APPROXIMATION 
1500  Meters 

Y^  = -9.7766  + .0198Y^  + (-.0218Yj^)  + (-.008Y^)  + 
.0033Yt^  + 1.4738Yj^ 

2500  Meters 

Y^  = -31.5553  + .0633Y^  + (-.1298Yjj)  + (-.0353Y^)  + 
.0153Yy  + 2.^553\ 

XM-1  FCS  APPROXIMATION 
1500  Meters 


Y^  = -9.7766  + .0198Y^  + -.0218(Yjj  ♦ .01)  + 

-.008(Y^  * .02)  + .0033(Y^  * .03)  + 1.4738Y^ 

2500  Meters 

Y^  = -31.5553+  .0633Yy  + -.1298(Yjj+  .01)  + 

-.0353(Y^  * .02)  + .0153(y,^  * .03)  + 2.4553Yj, 
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7.  Sununary 

This  section  has  presented  a basic  Miss  Distance  model  that  can 
be  used  to  determine  the  impact  point  of  each  round  fired  at  the  target 
plane,  either  in  single  shot  or  burst  mode  of  fire.  Hit  probabilities 
(P„)  have  not  been  computed  because  miss  distance  is  tlv,  basis  for  any 
later  generation  of  Pj^.  Although  the  model  as  presented  is  simplistic 
in  its  formulation,  since  much  of  the  laborious  computation  is  aleviated 
by  using  published  firing  table  data,  the  data  calculations  for  insertion 
into  the  model  variables  are  not  so  readily  available.  This  is  espe- 
cially true  when  attempting  to  numerically  define  tube  (barrel)  bend, 
effect  on  the  muzzle  due  to  trunnion  forces,  and  nonstandard  muzzle 
velocity  changes  due  to  nonstandard  propellant  temperature  and  barrel 
bend.  The  examination  of  all  of  these  input  variables  is  beyond  the 
scope  of  this  study,  however,  the  remaining  portion  of  this  chapter  will 
be  dedicated  to  an  in  depth  study  of  barrel  bend  and  its  resultant  im- 
pact demonstrated  in  Chapter  IV,  The  Simulation  Model. 


B.  BARREL  BEND 


1 . Introduction 

It  has  long  been  observed  that  if  a strip  of  metal  is  heated  and 
then  allowed  to  cool  in  such  a manner  that  one  side  cools  more  q.uickly 
than  the  other,  the  strip  will  bend  toward  the  cooler  side.  A more 
accurate  description  of  this  process  would  be  to  say  that  the  cooler 
side  returns  to  its  original  temperature  and  shape  more  q.uickly. 

This  process  also  occurs  in  gun  barrels.  When  the  gun  is  fired,  the 
barrel  is  heated.  The  barrel  is  exposed  to  the  environment  as  it  cools. 
If  sides  of  the  barrel  are  subjected  to  varying  environmental  conditions, 
e.g.  wind  or  rain,  the  barrel  will  cool  unevenly.  This  uneven  cooling 
causes  a temperature  differential  across  the  barrel.  The  magnitude  of 
this  temperature  differential  increases  over  time  to  a maximum  value 
and  then  diminishes  as  ambient  air  temperatirre  is  approached,  as  illus- 
trated in  Figure  3* 


Temperature 

T 

am 


Time 

T^  = Temperature  of  hotter  side 

Tj^  = Temperature  of  cooler  side 

T = Ambient  air  temperature 
am 

Figure  3s  Time-Temperature  Curves 

UO 


To  visualize  the  establishment  of  this  differential  more  clearly, 
imagine  that  the  temperature  of  all  points  of  a gun  barrel  cross  section 
has  been  raised  to  a pairticular  value . Suppose  further  that  a wind  is 
blowing  at  point  A in  the  figure  below. 


B 


D 


Figure  4:  Gun  Barrel  Gross  Section 

Point  A will  receive  the  full  benefit  of  the  cooling  effect  of 
the  wind  and  cool  much  more  qiiickly  than  Point  C.  Moving  from  Point  A 
toward  Points  B or  D the  temperature  will  rise  until,  in  the  vicinity 
of  those  points,  it  equals  that  of  Point  G.  The  temperature  would  then 
remain  relatively  constant  until  Point  G is  reached.  A temperature 
curve  for  the  surface  would  appear  as  shown  in  Figure  5" 

Rain  will  cause  a similar  differential.  In  this  case  the  hottest 
point  would  likely  be  near  Point  D and  the  coolest  near  Point  3.  The 
sun  can  ailso  affect  a temperature  differential  by  adding  temperature  to 
the  top  of  the  barrel. 

If  there  is  a temperature  differential  along  a gun  barrel,  the 
barrel  will  bend  in  the  direction  of  the  cooler  side . If  a round  is 
fired  through  a bent  gun  barrel  there  will  be  an  effect  on  the  intended 
trajectory  of  the  round.  In  order  to  approach  the  problem  of 
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# A ^ ‘ 


A B or  D G 

Surface  Points 


= Temperature  at  Point  G 
T^-AT  = Temperature  at  Point  A 

Figure  5;  Temperature  Gurve 

determining  the  effects  on  accuracy,  it  is  necessary  to  determine  how 
hot  the  outer  surface  of  the  barrel  gets,  how  much  of  a differential  can 
occur  across  the  barrel,  and  how  much  bend  will  result. 


The  highest  temperature  on  the  outside  surface  of  a gun  barrel 
(T^)  is  a complex  function  of  many  factors  including  barrel  thickness, 
material,  type  of  propellant,  number  of  rounds  fired,  rate  of  fire,  etc. 
T^  is  also  unique  for  each  different  type  of  gun  and  will  vary  slightly 
between  similar  guns.  Heat  and  temperatiure  curves  for  matny  guns  have 
been  developed.  Since  the  75^  cannon  has  been  selected  for  this 
analysis,  T^  must  be  approximated  for  this  particular  barrel. 
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AVERAGE  TEMPERATURE. 


RESIDUAL  BORE 
TEMPERATURE 


I /: 

; / 


OUTER  BARREL 

temperature 


rounds  fired 


6:  APE?  T'^mpera.turr'  Cur/OG 
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The  area  of  concern  is  near  the  muzzle  end  of  the  and  can  be 


assumed  to  be  representative  of  the  temperature  along  the  forward  end  of 
the  barrel.  The  outer  barrel  temperature  is  the  quantity  of  interest 
for  the  ainalysis.  The  curve  in  Figure  6 was  developed  by  a computer 
program  and  the  exact  mathematical  basis  could  not  be  provided.  The 
shape  of  the  curve  indicates  that  a linear  approximation  should  be  very 
satisfactory  for  estimating  T^  in  the  region  of  interest.  Linear  re- 
gression on  data  points  of  the  curve  results  in  the  following  approx- 
imation for  T . 

o 

Tq  = 64.55  + 19  (15) 


where  X = number  of  rounds  fired  at  maximum  rate . 

Figure  7 shows  the  data  from  the  Figure  6 curve  in  the  first 
two  columns.  T^  estimated  from  (15)  is  shown  in  the  third  column.  In 
the  last  column  is  T^  estimated  by  substituting  the  ambient  air  temper- 
ature (T  ) for  the  value  of  64.55  in  (15)* 

3.01  \ -"z 


Rounds  Regression  Regression 

Fired  Figure  6 (with  64,55)  (with  T.  ) 


0 

75 

64.55 

75 

5 

125 

I63.OO 

174 

10 

225 

262.00 

273 

15 

360 

362.00 

372 

20 

485 

460.00 

471 

25 

595 

559.00 

570 

30 

690 

658.00 

669 

35 

780 

757-00 

768 

40 

860 

856.00 

867 

45 

940 

956.00 

966 

50 

1010 

1055.00 

1065 

Figure  7:  Approximation  of  T^ 
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The  original  regression  and  the  regression  with  the  substitution 

of  T are  adequate  approximations,  especiailly  in  the  range  of  10  to  45 
3»in 

rounds.  This  approximation  is  used  to  estimate  the  temperature  differ- 
ential in  the  development  of  subsequent  relationships.  Trials  with 
differing  ambient  temperatures  show  that  while  T^  varies,  as  would  be 
expected,  the  temperature  differential  remains  the  same.  Also  the  use 
of  the  approximation  with  T^^  included  simplifies  some  later  calcu- 
lations.  The  coefficient  of  determination  (R  ) of  the  regression  in- 
cluding T^^  is  greater  than  .99.  Therefore,  for  the  balance  of  the  com- 
putations made,  T^  will  be  approximated  by  the  following: 

T = T + 19. 8X  (16) 

o am 

’(rtiere 

Tam  = ambient  air  temperature 

X = number  of  rounds  fired  at  maximum  rate 


3.  Approximating  the  Tenrperature  Differential 


A formula  for  the  effect  of  wind  or  rain  will  taJce  the  general 


AT  = G(T  - T ) 0<C  <1 

' 0 am"^ 


(17) 


where 

AT  = the  temperature  differential 
C = cooling  factor 

T^  = temperature  of  the  hottest  part  of  barrel  surface 
T = ambient  air  temperature 


45 


w 


In  other  words,  the  temperature  differential  that  will  develop 
across  the  barrel  will  be  some  fraction  of  the  difference  between  the 
hottest  point  of  the  barrel  and  the  ambient  air  temperature.  The  cool- 
ing factor,  C,  will  be  the  value  attainable  by  some  function  of  the 
intensity  of  wind  or  rain,  time  amd  other  factors.  The  search  for 
functions  representing  the  cooling  effect  of  wind  or  rain  was  most  un- 
successful. Neither  functions  nor  sufficient  data  to  permit  derivation 
could  be  found.  The  search  was  then  redirected  to  find  the  maxiraum  ex- 
pected value  for  C.  A value  of  G could  not  be  determined  for  rain. 
However,  determination  of  the  maximum  expected  value  of  G is  possible 
for  wind  cooling. 

By  Newton's  Law  of  Gooling: 

Q = h A 2iT  (18) 

where 

Q = heat  transfer  rate 

h = heat  transfer  coefficient 

at  = temperature  differential  between  the 
surface  and  ambient  air  temperature 

Ne  can  reasonably  assume  the  heat  transfer  rate  (Q)  on  the  sur- 
face of  a gun  barrel  to  be  constant.  The  definition  of  the  temperature 
differential  is: 

AT  = (T  - T^  J (19) 

3>iil 

We  cam  combine  (18)  and  (19)  under  the  above  assumption  to  give; 


Figure  8;  Gun  Barrel  Cross  Section 

Since  the  temperatiare  differential  between  T.  and  T in 

1 3>ni 

Figure  8 is  different  from  the  temperature  differential  between  and 
and  since  both  quantities  equal  a constant  value,  then  the  heat 
transfer  coefficient  (h)  must  vary  and: 

(T.  - T ) h = -^  (21) 

' 1 am  max  A ' 

(T  - T ) h , = (22) 

' o am^  min  A ' 

From  (21)  and  (22): 

(T,  - T ) h = (T  - T ) h . (23) 

1 am"^  max  o am*^  min 

Since,  by  definition: 

T.  = T - AT 
1 0 

(23)  can  be  reduced  to: 

- T 

ai 


(24) 


Solving  (25)  for  AT  yields* 

AT  " (1  - (T„  - Tj  (26) 

max 

The  ratio  of  h . and  h can  be  determined  from  the  work  of 
min  max 

Giedt.  [11]  With  a Reynolds  number  of  70,800  and  4.25  inch  diameter  gun 
barrel  the  value  of  the  ratio  is  found  to  be  IOO/25O  or  0.4.  Thus, 

the  wind  cooling  factor  for  the  ARES  cannon  is  0.6.  This  corresponds 
roughly  to  a 20  mph  wind.  Formula  (26)  for  the  temperature  differential 
may  now  be  reduced  to: 

where 


at  = the  maximum  expected  temperature  differential 
^ across  the  barrel  of  the  ARES  cannon  caused 
by  wind  cooling 

T^  = temperature  of  hottest  point  of  barrel  surface 
T = ambient  air  temperature 

With  the  above  expression  for  the  maocLmum  expected  temperature 
differential,  it  is  possible  to  calculate  the  magnitude  of  bend  that  can 
occur  in  the  ARES  cannon  ais  a result  of  wind  cooling.  This  requires  a 
formula  for  bend  which  is  developed  in  the  following  section.  Factually, 
this  is  the  only  bend  that  can  be  calculated  with  some  confidence.  It 
is  desirable,  however,  to  have  some  feel  for  the  range  of  bend  that  can 
occur  because  of  varying  conditions  of  rain,  as  well  as  wind.  There- 
fore, a parametric  analysis  was  performed  with  assumed  values  for  other 
conditions.  It  must  be  emphasized  that  the  value  of  0.6  for  wind  cool- 
ing is  the  only  value  that  has  been  derived.  This  value  will  be  used  in 
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the  study  as  the  cooling  factor  for  "hard  vrind."  Other  values  are  based 
solely  on  the  subjective  judgement  of  the  author  and  others  consulted.^ 


For  the  pvirpose  of  the  parametric  study  the  following  cooling  factor 

values  will  be  assumed  for  the  weather  conditions  indicatedi 

Hard  wind  - 0.6  (approximately  20  mph) 

Light  wind  - 0.2  (approximately  5 Kiph) 

Hard  rain  - 0.8  (approximately  0.5  inches  per  hour) 

Light  rain  - 0.6  (approximately  0.05  inches  per  hour) 

4.  Derivatj.on  of  the  Barrel  Bend  Formula 

From  the  definition  of  curvature: 


(28) 


where 

K = ciorvature 
9 = local  deviation  angle 
X = desired  axis 
it  can  be  seen  that 

0=^*  K dx  (29) 

The  task  is  to  find  an  expression  for  K.  Hooke's  Law  for  isotropic 
materials  I3  states; 

- Tj  (30) 

where 

£ ^ = strain  in  direction  of  length 


R.  E.  Newton  and  M.  D.  Kelleher  of  the  Department  of  Mechajiical 
Engineering,  Naval  Postgraduate  School  provided  assistance  here  and  in 
other  parts  of  the  derivation. 
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* A 


<T  = normal  stress  on  a barrel  cross  section 
X 

E = modulus  of  elasticity 

OL = coefficient  of  thermal  expansion 

T = temperature  of  a point 


T = average  temperat\ire 

ct 


Multiplying  (30)  by  y da  gives 

The  bending  moment  does  not  change  as  a result  of  temperature  strains 
cind  may  be  considered  zero  for  this  purpose.  Therefore, 

^^CTx  y da  = 0 (32) 

ajid  (31)  may  be  reduced  to; 


Since  by  Popov  C13] 


Ky 


(33)  becomes 


K Jy^  da  = fl^J(T  - T^)  y da 


By  definition; 


where 


y^  = -^  (r^^  - r^^) 


r = radius  of  outer  surface 
0 

r^  = radius  of  inner  surface 


* A . 


Substitution  of  (36)  into  (35)  yields: 


K = -f-  \ (T  - Tj  y da  (3?) 

It  is  necessary  to  maJce  certain  statements  and  assumptions  before 
proceeding  with  the  derivation.  If  a gun  barrel  cross  section  is  ex- 
periencing uneven  coolingi  a temperature  differential  will  develop.  As 
illustrated  in  Figure  9i  the  surface  temperature  at  0 = 0 will  be  equal, 
to  - AT.  The  temperature  will  rise  as  the  angle  0 increases  and 
reach  a maximum  value  near  0 = The  temperature  will  remain  con- 

stant xintil  0 = 7(.  Assume  the  point  of  equal  temperature  occurs  where 


Figure  9: 

Then  the  surface  temperature 


T 

T 


s 


s 


T 

0 


- AT  cos  0 


Cross  Section  Temperature 

(T  ) can  be  expressed  by  the  following; 
3 

0^O<-^  (38) 

01-^ 


(39) 


to  be; 


T = T - ^ In  (-|-)  (40) 

In  (— ) 
i 

Finally,  it  can  be  shown  that  relatively  little  bending  occurs 
at  the  breech  end  of  the  gun  baurrel.  This  fact  is  due  to  the  increased 
wall  thickness  and  because  of  the  various  supports  that  are  found  in  this 
portion  of  the  bairrel.  The  amount  of  bend  in  this  portion  is  assumed  to 
be  insignificant.  Therefore,  the  portion  of  the  barrel  to  be  considered 
is  the  portion  which  extends  forward  of  the  supported  section.  This 
portion  will  be  referred  to  as  the  unsupported  length.  It  is  also  noted 
that  in  the  unsupported  length  of  the  barrel,  the  inner  and  outer  radii 
are  generally  constant  for  most  guns,  including  the  ARES  cannon.  If  the 
radii  are  not  constant  in  this  region,  then  an  expression  must  be  in- 
serted in  the  relationships  which  define  the  radius  as  a function  of 
length.  Inner  and  outer  radii  will  be  considered  constant  through  the 
unsupported  length  of  the  bacrrel  for  the  remainder  of  the  derivation. 

From  (37)  let: 

P = J (T  - T^)  y da  (41) 

in  cylindrical  coordinates 


Substituting  (45)  and  (46)  into  (44)  gives 


Y cos  0 d0  -X\  (T.  - T ) cos  0 d0  (4?) 

0 Jo  ^ ® 


From  the  definitions  of  T given  by  (38)  and  (39).  (4?)  may  be  expanded 

s 


to: 


Integration  with  respect  to  0 gives 


p = 


(^9) 


(49)  can  now  be  replaced  in  (37)  to  give 

which  by  (36)  can  be  further  reduced  to 


A Of  At 
K - —5 JT 

r - r. 
o 1 

(51)  can  now  be  substituted  into  (29)  to  give; 

X d At 


0 


■i: 


IT 

r - r. 
o i 


dx 


(50) 


(51) 


(52) 


The  result  of  integration  of  (52)  will  be  in  radians.  Since  it 
is  desired  to  have  the  result  in  mils,  the  right  side  of  (52)  is  multi- 
plied by  the  conversion  factor  of  1018.6  mils  per  radian.  Since  the 
terms  in  (52)  are  constants,  including  the  mil  conversion  factor,  let 
a bending  constant  ( p ) be  defined  where: 


g _ (lQ18.6)(X)(d) 
r Zj; 


(53) 


Then  (52)  may  be  reduced  to 

0 = pAT  dx  (54) 

Since  the  terms  to  be  integrated  are  constant  along  the  unsupported 
length  of  the  barrel,  the  result  of  integration  will  give  the  amovmt  of 


bend  at  the  muzzle  in  mils  and  will  be; 
0 = L p AT 


(55) 


where ; 

L = length  of  unsupported  barrel  in  millimeters 
P = bending  constant  defined  by  (53) 

AT  = temperature  differential  defined  by  (26) 

5.  Calculating  Bend  for  the  ARES  Cannon 

Relationships  derived  in  the  foregoing  sections  permit  calcula- 
tion of  the  maximum  expected  bend  for  the  ARES  gun  as  a result  of  wind 
cooling.  For  the  ARES  cannon  [l^  s 
r^  = 5^  mm 
r^  = 37.5  nun 
L = 1700  ram 
= 6 X 10”^ 

Inserting  the  values  above  into  {^6)  , (53)  and  (55)  gives  the  final 
result: 


e = 0.0327  AT 


(56) 


(27)  shows  that  for  wind  cooling: 


AT  = 0.6(T  - T ) 

' o am' 


(57) 


By  (16),  T^  may  be  approximated  by: 


T = T + I9.8X 
o am 


(58) 


Substituting  (58)  into  (57)  and  solving  for  AT,  and  then  substituting 
the  result  into  (56)  yields 


e = 0.39X 


(59) 


where  X is  the  number  of  rounds  fired  at  maximum  rate. 
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1. 

[ 

This  leads  to  the  conclusion  that  the  maximum  expected  bend  at 
the  muzzle  of  the  ARES  cannon  as  a result  of  wind  cooling  can  be  approx- 
imated by  (59) . Recall  that  (58)  is  a good  approximation  especially 
in  the  range  of  10  to  45  rounds.  A graphic  representation  of  (59)  is 
given  in  Figure  10 . 

20 
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(mils) 
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Figure  10;  Wind  Generated  Bend  (ARES  Gun) 
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As  can  be  seen  from  Figure  10,  the  amoimt  of  bend  that  can  be 
caused  by  wind  cooling  can  be  very  significant.  These  results  are 
markedly  less  severe  than  has  been  estimated  by  other  researchers.  Rain 
can  cause  even  more  severe  bend. 

It  must  be  recognized  that  it  is  extremely  unlikely  that  a tank 
crew  would  fire  45  rounds  at  maximum  rate.  However,  given  the  capabil- 
ity, it  is  not  unlikely  that  10  rounds  would  be  fired.  Several  shorter 
bursts  fired  in  rapid  succession  will  generate  nearly  the  same  amount 
of  temperature  and  the  resultant  bend.  An  additional  question  must  now 
be  addressed.  Given  that  the  barrel  is  bent,  how  much  deflection  to  the 
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6 . Straightening  Forces 

If  a round  is  fired  through  a bent  barrel,  the  projectile  may 
exert  a straightening  force  on  the  barrel  as  it  travels  down  the  bore. 

If  there  is  a degree  of  straightening,  the  amoimt  of  actual  deflection 
will  be  less  than  the  full  amount  of  bend.  Any  amount  of  straightening 
that  may  occur  must  be  paid  for  by  an  accompanying  loss  in  muzzle  veloc- 
ity. A loss  in  muzzle  velocity  will  cause  a range  error.  The  round  will 
not  travel  as  far  as  intended  and  will  pass  through  the  target  plane  at 
a point  lower  than  planned.  There  is  no  easy  way  to  determine  the  net 
deflection  auid  any  chajige  in  muzzle  velocity. 

The  results  of  the  parametric  analysis  conducted  in  Chapter  IV 
indicate  that  the  amount  of  bend  will  seldom  exceed  one  degree.  This 
amount  of  bend  is  very  significant  as  far  as  accuracy  is  concerned,  but 
relatively  insignificant  with  respect  to  the  forces  involved.  It  is 
assumed  that  for  the  small  angles  encountered,  the  amount  of  straignten- 
ing  will  be  negligible.  The  full  amount  of  bend  will  then  be  considered 
as  deflection  of  the  round  from  its  intended  trajectory  and  there  will 
be  no  resulting  loss  of  muzzle  velocity. 

7 . Time  Dependency 

Figure  3 indicates  the  magnitude  of  bend  varied  over  time . The 
angle  will  increase  to  a maximum  magnitude  and  then  diminish  as  each 
side  of  the  barrel  approaches  the  ambient  air  temperature.  The  maximim 
expected  bend  from  wind  cooling  will  occur  approximately  15  minutes 
after  completion  of  firing.  A hard  rain  will  cause  the  maximum  bend  to 
occur  in  a much  shorter  time.  This  time  cannot  be  accurately  calculated 
due  to  the  lack  of  adequate  data.  It  is  assumed  in  this  study  that  the 
maximum  maignitude  of  rain  induced  bend  will  occur  approximately  five 
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minutes  after  completion  of  firing.  The  time  required  for  the  maudmum 
value  of  bend  to  occur  will  vaLry  between  types  of  guns  and  with  varying 
amounts  of  wind  or  rain. 


8 . Avoiding  Bairrel  Bend 

It  is  a fact  that  gun  barrels  bend.  The  foregoing  sections  have 
demonstrated  that  the  amount  of  bend  can  be  very  significamt  and  can 
cause  a sizable  miss  distance  at  the  target  plane.  This  effect  can  be 
avoided  by  two  general  methods.  The  amount  of  bend  occuring  at  any 
point  in  time  can  be  compensated  for  or  barrel  bend  can  be  eliminated 
altogether . 

Methods  of  compensating  for  barrel  bend  include; 

a.  Temperature  sensing  of  points  of  the  barrel  to  determine 
temperature  differential,  calculating  the  amount  of  bend,  and  adjusting 
the  lay  of  the  gun  accordingly. 

b.  Direct  measurement  of  bend  thro\igh  optic  or  other  sensing 
devices  and  compensating  as  above . 

There  are  several  ways  of  avoiding  barrel  bend  altogether. 

These  basically  involve  the  prevention  of  uneven  cooling.  Two  methods 


a.  Installation  of  a cooling  system  ;diich  insures  even  cooling 
of  all  parts  of  the  barrel. 

b.  Shield  the  barrel  from  the  causes  of  uneven  cooling  (i.e.  the 
environment)  with  a protective  shroud  or  coating.  This  alternative  is 
reportedly  being  developed  by  certain  allied  nations. 

9.  Summary 

Gun  barrels  are  heated  dvaring  firing.  Rapid  fire  guns  in  partic- 
ular can  generate  rather  high  barrel  temperature,  although  high  barrel 
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siirface  temperatures  are  certainly  common  in  guns  vfith  slower  rates  of 
fire  when  a large  number  of  rounds  are  fired.  The  gun  barrel  is  ex- 
posed to  varying  environmental  effects  as  it  cools  after  completion  of 
firing.  Wind  or  rain  on  the  surface  of  the  barrel  during  the  cooling 
period  will  cause  the  barrel  to  cool  unevenly,  resulting  in  a bend  of 
the  barrel.  This  occurrence  can  significantly  affect  the  trajectory  of 
rounds  fired  while  the  bairrel  is  bent  causing  a sizable  miss  distance 
at  the  plane  of  the  target.  The  effects  of  barrel  bend  can  be  compen- 
sated for  or  barrel  bend  can  be  avoided  altogether.  Weapons  under 
development  should  be  required  to  incorporate  methods  of  avoiding  barrel 
bend  or  its  effects  and  modifications  to  existing  weapons  should  be 
investigated. 


t 
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IV.  THE  SIMULATION  MODEL 


A.  MODEL  DESCRIPTION  AND  ASSUMPTIONS 

To  employ  the  miss  distance  and  barrel  bend  models  under  Quasi- 
Combat  conditions,  a computerized  parametric  analysis  was  conducted.  It 
was  performed  in  order  to  evaluate  the  automatic  gun  concept  under 
various  weather  conditions,  fire  control  system  configurations  and  pro- 
jectile types.  When  required,  the  analysis  used  specific  characteristics 
of  the  ARES  75^  automatic  gun. 

To  evaluate  the  effects  of  weather  on  gun  performance  (and  in  par- 
ticular the  effects  of  barrel  bend)  a weather  case  matrix  was  developed 
containing  variable  inputs  of  air  temperature,  direct  left  to  right 
crosswind,  rainfall,  vapor  pressure,  and  cooling  factors  which  were 
developed  in  Chapter  III,  Section  B.  The  weather  condition  for  any 
trial  is  considered  to  occur  at  1000  feet  above  sea  level.  Air  temper- 
ature varies  from  the  standard  of  60°  Fahrenheit  to  either  100°F  or  0°F. 
Rangewind  has  not  been  simulated  but  direct  left  to  right  crosswind 
effects  can  be  determined  by  using  a light  wind  of  approximately  7 inph 
and  a haxd  wind  of  approximately  20  mph.  The  effects  of  rain  are  simu- 
lated by  considering  a light  rain  condition  of  0.05  inches  per  hour  and 
a hard  rain  condition  of  0.5  inches  per  hour.  Vapor  pressure  values 
(in  inches  of  mercury)  aure  a function  of  the  air  temperature  and  relative 
humidity  and  are  used  in  the  model  to  compute  air  density.  For  the  "no 
rain"  condition,  the  relative  humidity  was  assumed  to  be  80^.  If  a rain 
condition  existed,  then  the  relative  humidity  was  set  at  100^.  By 


referring  to  vapor  pressure  tables  in  Smithsonian  Institute  ^Sj,  spe- 
cific values  were  then  generated  for  vapor  pressure  considering  the  air 
tempera tvire  and  relative  humidity. 

Each  test  weather  case  was  conducted  in  conjunction  with  a selected 
projectile  (HE  or  APFSDS)  and  range  to  engage  the  target.  Once  all 
selected  weather  cases  had  been  simulated  for  a specific  test  projectile 
and  range,  then  all  weather  cases  were  repeated  with  the  next  test  pro- 
jectile and  range . 

After  generating  the  basic  situation  by  establishing  a weather  case, 
the  test  projectile  and  range  of  the  target,  the  only  element  remaining 
in  the  model  scenario  is  the  trial  itself.  The  basis  for  the  structure 
of  a trial  is  an  attempt  to  structure  a somewhat  realistic  combat 
scenario  for  a rapid  fire  gun  mounted  on  a tank. 

The  trial  represents  perhaps  an  hour  in  combat  for  a tank.  It  is 
envisioned  that  during  this  period  a tank  would  have  several  short 
engagements . During  each  engagement  the  tank  would  fire  several  short 
bursts  of  either  three  or  five  rounds  at  the  target.  Between  engage- 
ments the  tank  could  remain  in  position  or  move  to  a new  location.  For 
this  study  three  engagements  per  trial  and  four  bursts  per  engagement 
were  used.  (See  Figure  11) 

The  time  between  engagements  is  established  as  the  time  required  for 
the  maximum  expected  bend  to  occur.  This  time  will  vary  according  to 
the  environmental  conditions  being  experienced  by  the  barrel.  Maximum 
bend  can  occur  in  as  little  as  five  minutes  under  a hard  rain  condition 
or  longer  for  lessor  weather  conditions.  The  length  of  any  engagement 
is  less  than  one  minute . An  engagement  is  thus  considered  a sequence 
of  four  short  bursts  fired  in  rapid  succession  at  maximum  rate  of  fire 
with  minimal  time  between  bursts. 
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ENGAGEMENT 


BURST 


ROUNDS/BURST 


3 or  5 

3 or  5 
3 or  5 
3 or  5 

3 or  5 
3 or  5 
3 or  5 
3 or  5 
3 or  5 
3 or  5 

3 or  5 
3 or  5 


Figure  111  Trial  Representation 

The  gun  barrel  surface  heat  generated  in  an  engaigement  will  nearly 
approach  that  generated  by  firing  a single  burst  of  the  same  number  of 
rounds.  The  program  structxare  calls  for  the  gun  barrel  surface  to  be 
at  ambient  air  temperature  for  initialization  of  each  trial.  Therefore, 
there  is  no  bairrel  bend  for  the  first  engagement.  The  number  of  rounds 
fired  during  the  first  engagement  is  used  to  calculate  the  resulting 
gun  barrel  surface  temperature.  Bend  at  the  muzzle  in  mils  is  then 
calculated  depending  on  the  weather  case  selected.  Since  wind  is 
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assumed  to  be  blowing  left  to  right  and  perpendicular  to  the  barrel,  | 

I 

I 

bend  in  the  X aois  is  always  to  the  left.  Rain  is  assumed  to  be  falling  | 

on  top  of  the  bsmrel.  Rain  induced  bend  is  calculated  in  exactly  the  j 

i 

same  manner  with  only  the  cooling  factor  changed.  Thus,  bend  in  the  Y 
direction  is  always  upward. 

For  the  third  engagement  of  a trial  one  final  factor  must  be  con- 
sidered. The  number  of  rounds  in  the  second  engagement  generate  barrel 
surface  temperatiire  which  is  added  to  the  temperature  remaining  from  the 
first  engagement.  Thus  the  bend  generated  between  the  second  and  third 
engagements  will  be  more  severe . After  consideration  of  the  time  span 
and  the  foregoing  assumptions,  it  was  determined  that  a standard  cumu- 
lative heat  value  of  0.5  should  be  used.  Therefore,  the  baumrel  surface 
temperature  at  the  end  of  the  second  engagement  is  based  on  the  number 
of  rounds  in  the  second  engagement  plus  O.5  times  the  number  of  rounds 
in  the  first  engagement. 

B.  MODEL  RESULTS 

Output  tables  from  the  computer  trials  are  given  in  Appendix  B.  In 
order  to  simulate  the  two  systems,  the  standard  deviations  developed  in 
Chapter  III,  Section  5 >fere  used.  The  basic  differences  in  the  two 
system  approximations  were  due  primarily  to  the  higher  sophistication  of 
the  XM-1  fire  control  employing  sensor  devices  and  the  laser  rangefinder. 

Table  VII  is  a sample  of  the  computer  output  to  demonstrate  the 
model's  capabilities.  One  trial  is  presented  and  the  applicable  general 
trial  situation  is  listed  in  the  page  heading.  In  this  case,  the  trial 
situation  occurs  using  the  ARPA  HE  round  and  engsiging  a target  at  a 
range  of  2500  meters  under  weather  conditions  consisting  of  an  air 
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1 

temperature  of  0°F  vdth  a hard  orosswind  of  approximately  20  mph  and  no 
rain.  For  reader  clarity,  some  of  the  column  headings  presented  for  the 
trial  results  are  defined  helow. 

1 . FIXED  BIAS  (mils)  is  the  deviation  in  mils  at  the  muzzle  of  the 
fixed  biases  of  the  weapon  system  for  any  round  fired. 

2.  BEND  is  the  amount  of  barrel  bend  on  the  muzzle  expressed  in  mils 
under  the  scenario  conditions. 

3.  TOTAL  is  the  total  fixed  bias  in  mils  at  the  muzzle  under  the 
scenario  conditions,  including  barrel  bend. 

4.  RANDOM  ERROR  (MILS)  is  the  effect  of  Round- to-Round,  Burst- to-Burst, 

and  Occasion- to-Occasion  errors  on  the  projectile  expressed  in  mils  at 
the  muzzle.  These  error  values  were  generated  using  the  error  classi- 
fications and  computation  of  standard  deviations  as  presented  in  Chapter 
III  and  then  applying  the  error  values  by  use  of  a standard  normal  ran- 
dom number  generator  routine.  After  completion  of  the  computer  runs,  a j 

potential  problem  with  the  random  number  routine  was  discovered.  There-  I 

fore,  this  routine  should  be  completely  verified  before  using  the  model  ‘ 

in  a production  mode . 

5.  AZIMUTH  refers  to  the  X axis  of  the  target  plane. 

6 . EIEVATION  refers  to  the  Y axis  of  the  target  plane . 

7.  TOTAL  ERROR  (METERS)  is  the  miss  distance  in  meters  of  the  projec- 

I 

tile  at  the  target  due  to  the  contribution  of  all  the  fixed  biases, 
variable  biases  and  random  errors. 

As  presented  in  the  body  of  the  table,  the  trial  consists  of  three  i 

) 

engagements,  each  less  than  one  minute  in  duration,  foirr  short  bursts  i 

fired  in  rapid  succession,  and  three  rounds  per  burst.  For  Engagement  1, 
gun  barrel  surface  temperature  is  at  ambient  air  temperature  and 

:‘j 
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therefore,  barrel  bend  equals  0.0.  The  factor  affecting  azimuth  fixed 
bias  in  this  situation  is  the  hard  left  to  right  crosswind  resulting  in 
3.07  mils  bias  at  the  muzzle.  Factors  which  affect  the  elevation  fixed 
bias  are  the  temperatvire  differential  from  standard  of  60°F  to  0°F  and 
the  air  density  change  from  standard  due  to  the  cold  temperature . The 
random  error  totals  at  the  muzzle  are  the  sum  of  the  Round- to-Round, 
Burst- to-Burst,  and  Occasion- to-Occasion  errors  after  use  of  the  standard 
normal  random  number  generator  routine.  Notice  that  the  Occasion-to- 
Occasion  errors  do  not  vary  within  the  engagement;  the  Burst- to-Burst 
errors  vary  only  between  bursts  within  the  engagement;  and  the  Round-to- 
Round  errors  vary  with  each  round  fired.  Finally,  the  total  miss  dis- 
tance in  meters  at  the  target  is  expressed  for  each  round  in  azimuth  and 
elevation  by  considering  both  the  system  fixed  bias  and  the  random 
errors . 

In  Engagement  2,  the  factors  affecting  the  fixed  biases  of  the  wea- 
pon system  remain  the  same;  however,  now  barrel  bend  is  introduced  as  an 
additional  fixed  bias.  The  12  rounds  fired  in  the  first  engagement 
affect  the  gun  barrel  surface  temperature  and  the  time  between  the  two 
engagements,  defined  as  the  time  for  maximum  bend  to  occur,  has  allowed 
the  cool  hard  wind  to  decrease  the  outer  surface  temperature  with  the 
resultant  bend.  The  new  fixed  bias  in  azimuth  becomes  4.6?^  and  7.74^ 
for  barrel  bend  and  total  fixed  bias  respectively. 

In  Engagement  3i  barrel  bend  is  the  only  fixed  bias  change.  The 
number  of  rounds  fired  in  the  first  engagement  (12)  and  the  rounds 
fired  in  the  second  engagement  (12)  are  cumulatively  ad.dressed  to  deter- 
mine the  resulting  outer  surface  temperature  of  the  barrel.  Although  the 
same  crosswind  velocity  is  encountered,  the  effect  of  the  higher  outer 
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surface  barrel  temperature  is  a more  severe  barrel  bend  condition  of 
7.01/.  This  results  in  a total  fixed  bias  in  azimuth  of  10.08;il. 


There  were  a total  of  480  trial  runs  conducted  in  the  simulation 
model  (240  trials  with  M6O  fire  control  system  approximation  and  240 
trials  with  XM-1  fire  control  system  approximation) . The  results  of  the 
selected  trial  runs  included  in  Appendix  B adequately  illustrate  the 
effects  of  barrel  bend  in  automatic  cannon  firing  and  validate  the 
models  developed.  If  additional  trial  run  output  is  required  for 
analysis,  they  are  available  upon  request  from  the  authors. 
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SIMMAHY  AND  CONCLUSIONS 


Conventional  warfare  is  a realistic  possibility  in  the  turbulent 
world  of  today.  In  a conventional  war  the  tank  will  be  a principle 
threat  to  allied  forces  and  the  primary  counter  to  that  threat.  If  the 
tank  is  to  be  effective  in  the  antitank  role,  in  view  of  the  numerical 
superiority  of  opposing  forces,  there  must  be  a capability  of  accurate 
fire  well  above  that  demonstrated  in  the  Arab-Israeli  War.  Every  attempt 
must  be  made  to  eliminate  or  minimize  factors  that  affect  the  accuracy 
of  tank  gun  fire . 

This  paper  has  presented  an  analysis  of  factors  affecting  miss 
distance  and  proposes  that  barrel  bend  has  been  overlooked  ais  a major 
contributor  to  inaccurate  fire . A formula  was  derived  for  computing 
the  expected  value  of  bend  which  can  realistically  be  experienced  in 
varying  weather  conditions.  A model  was  presented  which  computes  barrel 
bend  and  other  errors  giving  the  resulting  miss  distance  at  the  target. 

Analysis  of  the  model  output  as  presented  in  Appendix  B provides 
the  following  conclusions. 

1.  Miss  distance  at  the  target  is  less  with  the  XM-1  fire  control  system 
approximation  than  the  M60  fire  control  system  approximation  due  primar- 
ily to  the  capability  of  sensors  for  air  temperature,  wind,  and  air 
density  as  well  as  the  laser  rangefinder.  This  minimizes  significantly 
the  system  fixed  biases  other  than  barrel  bend.  The  m60  fire  control 
system  is  particularly  susceptible  to  error  induced  by  crosswinds. 

The  fixed  bias  unaccounted  for  by  the  XM-1  fire  control  was  generally 
less  than  0.1  mil.  This  error  is  relatively  insignificant  when  compared 
to  the  error  from  barrel  bend. 
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2.  The  APFSDS  round  is  more  accurate  than  the  HE  round  as  a result  of 
the  projectile's  design  characteristics,  e.g.  significantly  higher  muzzle 
velocity. 

3.  Rapid  fire  cannons  will  he  particularly  susceptible  to  barrel  bend 
due  to  the  high  barrel  stirface  temperatures  resulting  from  biirst  fire. 

A-.  Shorter  bursts  generate  less  heat  and  less  barrel  bend. 

5.  Cooler  air  temperatures  increase  the  magnitude  of  barrel  bend. 

6.  Shorter  barrels  and  barrels  with  thicker  walls  will  be  less  sus- 
ceptible to  bend. 

7.  Barrel  bend  occurs  during  the  cooling  period  after  an  engagement. 

The  magnitude  of  bend  that  can  occur  depends  on  the  number  of  rounds 
fired  during  the  engagement,  weather,  and  gun  characteristics.  An  en- 
gagement of  20  rounds,  i.e.  four  bursts  of  five  rounds  fired  in  rapid 
succession,  is  shown  to  generate  sufficient  barrel  temperature  to  cause 
as  much  as  a 12  mil  error  under  wind  cooling  or  nearly  16  mils  under 
rain  cooling.  This  bend  becomes  a fixed  bias  affecting  the  accuracy 

of  rounds  fired  in  the  next  engagement. 

8.  In  a series  of  engagements  there  will  be  an  additive  effect  of 
barrel  bend  if  sufficient  time  is  not  allowed  for  the  bend  to  diminish. 
Thus  the  third  engcLgement  can  be  affected  by  the  bend  induced  by  the 
second  engagement  plus  a portion  of  bend  induced  by  the  first  engaigement. 

There  are  ways  to  compensate  for  barrel  bend  and  ways  to  eliminate 
bend  altogether.  If  compensation  for  barrel  bend  is  to  be  achieved, 
then  extensive  research  is  required  toj 

1 . Develop  functions  for  outside  barrel  temperature . 

2.  Develop  functions  for  the  cooling  factor  of  varying  amounts  of  wind 

ckiiu  xalii  • 
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3.  Determine  whether  the  barrel  does  or  does  not  straighten  by  any 
significauat  amount  during  firing,  and  if  straightening  does  occur,  at 
what  cost  in  muzzle  velocity  from  standard. 

4.  Develop  systems  for  sensing  or  direct  measurement  of  barrel  bend. 

On  the  other  hand,  if  elimination  of  bend  is  desired,  then  research  is 
required  to  develop  a system  which  will  insure  even  cooling.  Shielding 
of  the  barrel  from  the  elements  of  wind  and  rain  by  a protective  shroud 
or  coating  seems  to  be  the  most  viable  alternative.  This  method  is  re- 
portedly being  investigated  by  West  Germany  and  Great  Britain. 

In  addition  to  further  analysis  of  barrel  bend,  the  testing  of  the 
models  in  the  parametric  analysis  presented  further  areas  for  extended 
study . 

1.  Examination  of  the  effects  of  trunnion  forces  on  muzzle  position  to 
provide  additional  input  into  the  miss  distance  model. 

2.  Determine  the  change  in  muzzle  velocity  from  standard  as  a result 

of  changes  in  propellant  temperature  due  to  heat  transfer  from  the 

barrel  temperature. 

3.  Determine  the  change  in  muzzle  velocity  from  standaird  as  a result 

of  the  barrel  bend  condition  and  projectile  straightening  upon  firing. 

This  paper  has  particularly  addressed  the  bend  that  can  occur  in 
rapid  fire  cannons  because  of  the  higher  barrel  surface  temperatures 
that  are  likely  to  be  experienced  as  these  weapons  are  developed.  The 
paper  also  concentrates  on  tank  guns  due  to  the  stringent  requirement 
for  accuracy  when  firing  at  point  targets  with  direct  fire.  However, 
bend  occurs  in  weapons  with  lower  rates  of  fire,  including  current  tank 
guns  and  conventional  artillery  weapons. 

It  is  concluded  that  barrel  bend  can  significantly  decrease  the 
accuracy  of  gunfire.  The  magnitude  of  the  resultant  error  at  the  target 

70 


I 


4 


can  no  longer  te  tolerated  if  the  tank  is  to  remain  the  principal 
anti-tank  weapon  on  the  modem  battlefield.  Weapons  under  development, 
particularly  rapid  fire  tank  guns,  should  be  required  to  demonstrate 
the  capability  of  avoiding  or  at  least  compensating  for  barrel  bend. 


APPENDIX  A 


ASSUMED  FIRING  TABLE  AND  EFFECTS  DATA 

The  firing  table  data  presented  in  Appendix  A is  not  the  exact  data 
being  used  for  the  ARPA  75”™  cannon  and  selected  ammunition,  but  is 
merely  aui  assumed  approximation  for  use  in  the  models  generated  in  this 
paper.  For  more  effective  evaluation  of  this  developmental  system,  the 
exact  firing  table  data  should  be  used. 


ASSUMED  FIRING  TABLE  AND  EFFECTS  DATA 
ARPA  APFSDS  - 75MM 


RANGE 

SUPEREIEVATION 

DRIFT 

(Meters) 

(MILS) 

(MILS) 

1500 

2.8370 

0 

2500 

4.8929 

0 

NON  STANDARD  - STANDARD 

RANGE 

CHANGE  IN 

CHANGE  IN 

NON-STD 

JuL 

HEIGHT  (M) 

DRIFT  (m) 

CONDITION 

1500 

.0060 

.0000 

1 m/s  Velocity 

-.0025 

.0000 

1 Pet  Density 

-.0015 

.0000 

1 Pet  Air  Temp 

.0000 

.0000 

1 m/s  Rangewind 

.0000 

.0358 

1 m/s  Crosswind 

1 .4730 

.0000 

1 MIL  Elevation 

2500 

.0144 

.0000 

1 m/s  Velocity 

-.0110 

.0000 

1 Pet  Density 

-.0030 

.0000 

1 Pet  Air  Temp 

.0015 

.0000 

1 m/s  Rangewind 

.0000 

.0983 

1 nys  Crosswind 

2.4540 

.0000 

1 MIL  Elevation 

c 


ASSUMED  FIRING  TABLE  AND  EFFECTS  DATA 
ARPA  HE  AMMO  - 75MM 

DRIFT 
(MILS) 

0 

0 

NON  STANDARD  - STANDARD 


RANGE 

(Meters) 

1500 


2500 


SUPERELEVATION 

(MILS) 

6.6336 

12.8519 


RANGE 

CHANGE  IN 

CHANGE  IN 

NON-STD 

(M). 

HEIGHT  (m) 

DRIFT  (m) 

CONDITICN 

1500 

.0198 

.0000 

1 m/s  Velocity 

-.0218 

.0000 

1 Pet  Density 

-.0080 

.0000 

1 Pet  Air  Temp 

.0033 

.0000 

1 m/s  Rangewind 

.0000 

.1970 

1 ra/s  Crosswind 

1.4738 

.0000 

1 MIL  Elevation 

2500 

.0633 

.0000 

1 m/s  Velocity 

-.1298 

.0000 

1 Pet  Density 

-.0353 

.0000 

1 Pet  Air  Temp 

.0153 

.0000 

1 m/s  Range  vrind 

.0000 

.6140 

1 m/s  Crosswind 

2.4553 

.0000 

1 MIL  Elevation 
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APPENDIX  B 


OUTPUT  TABLES 

The  computer  trial  output  tables  included  in  this  appendix  represent 
a selected  sample  of  the  model  output  for  analysis.  The  effect  of  each 
of  the  four  selected  weather  conditions,  vrtiich  are  constant  for  each  of 
the  eight  tables,  are  examined  by  varying  the  trial  conditions  of  the 
round,  range,  fire  control  and  rounds  in  the  burst. 
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HAi<D  RAIN  5 0.00  I 0.00  I 0.00  1 0.08  I 0.00  i 0.00  10.38  |10  .46  I 0.00 


HARD  RAIN  5 0.00  0.00  0.00  0.34  0.00  I 0.00  |l0.38  10.?2  0.00 


HARD  RAIN  5 0.00  0.00  0.00  0.00  O.OOl  0.00  1 10.38  1 0. 38  0.00 
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APPENDIX  C 


COMPUTER  PROGRAM 

The  computer  program  included  in  this  appendix  is  that  written  for 
simulation  of  the  miss  distance  and  barrel  bend  models.  The  only  changes 
required  to  simulate  the  various  fire  control  system  approximations 
(M60  or  XM-l)  are  to  change  the  data  inputs  for  percentage  of  residual 
sensing  errors  and  for  Round -to -Round,  3urst-to-Burst,  and  Occasion- to- 
Occastion  dispersion.  The  program  may  also  be  adapted  to  different  gun 
barrels  by  inserting  the  physical  characteristics  of  the  g\in  and  an 
approximation  for  gun  barrel  surface  temperature. 
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iciitnr>ooonov'T|nofhr»r»ooiioon>r»n^onior>rhooni'^n(or>(^r»r^ 


iCTITY-CALrULATEn  A I HI^NSITY 
-AC  IB-WC^ACXUM— l-NC- 


AC''VEl=CALCULAT'^n  muZZFL  VELOCITY 

ALF^»LENGTH  OF  UNSUPPCBT EO  onRTICN  OF  PARREL  IN  YM 

AtRMA=CO?FFICteNT  OF  THERMAL  EYPANSICA  

ALTOEF-EIRING  table  effect  hf  I MIL  MIJZZL=  LATERAL  OEFL'CTION 

ATMP=AIR  TFMPFRATijpf 

'IC^^=Cr^STA^T  SEND  COFFFICI=f)r 

RTeX=eURST  TO  BURST  AZIMUTH  STANOARC  CEVIATICN 
-aie-Y^SLiRST  TfT  BURST  F vFvAJ-^X^^-S-IA^40ARC-GeV-^ A-I-tgA 


CCM  = AIR  OFNSITY  RESinUAL  SENSING  EHRCR 

ccaz’Air  tempfratupe  resioual  sensing  error 
'CeNJ=PiND  resioual  SEPSING  =RR0R 

r.UPHT  = CUMULATIvE  pFAT  FACTOR 

CVihC=CRCSSW  IMO  left  to  RIGHT 
OFLTX^TEHPF BATURE  0TFFER=NTIAL  IN 
oeltystfmpfraturf  n I f fer fnt I al  in 
-Ct-CCCll  N C-  F ACT-GR- 


AXIS 

AXIS 


C 

r 

C 

c 

c 

-C 


FTATmRsFIriag  table  effect  GF  IX  CHANGE  FROM  STANOAPO  AIR  TEMP 
ETCWNC*FIRI AG  table  EFFECT  OF  C OOS S W INC ( M/ S I 
— F TONTY^FIRI AG  TARL=  EFFECT  OF  I*  CHANGE  FROM  STANCARO  AIR  CENSITY 
FirRFT^FIPINC  TABL=  PRCJECTILE  DRIFT 

FTFLfv=FIRIAG  TABL^  effect  of  I MIL  ELEVATION  OF  “UZZLF 
TTRCWDSFtRIAG'TASLS  EEFECr  gif  B ANGE  ti  I A C ( m/ S ) - 

E tfnge=fip.  1 NO  table  range 

E_L4PS4^;r.I  NO— lABCs— SUP!GaFWeVJAJLlC.N 

FTTPROsF IR I NG  table  TYPE  ROUND 

ETVLTY=FIRING  TABL=  fFFFCT  of  muzzle  VELOCITY  CFANGE  FROM  STAACARC 

MPESTiAL-MBFr  of-bursts  

nfag=fagagement  number 

ARCS'AUMBER  of  PCUNOS 

— mrxsnccssjcN  to  occasion  azimutp" staacard  devuticn 

rTCY=OCCASI ON  TO  OCCASION  EIFVATION  STANOARO  CEVIATION 

B a lAaBAlA  

bI=IANER  barrel  RADIUS  THROUGH  UNSUPPCRTEO  LENGTH  IN  MM 
oc  = cutfp  barrel  raoius  ’’hrough  unsupported  LENGIH  in  mm 

RTPX^RCUNO  TO  RCUNO  AZIMUTH  STANCARC  CEVIATICN  - - 

RTPY=RrUNO  TO  RCUMO  EL-VATION  STANOARC  OFVIATICA 

statmp^stancarc  air  tempfrature  at  1000  feet  _ 

'TCHTy^STANOARC^ A IR  OEASf TV  AT  1000  FEET  ’ 

5TCVEL=EIRING  table  STO  muzzle  VELOCITY 

\i.sracc-\^Aonp op  c c e i fp  = 

VflBB’CHANGE  in  ^LOCITY  from  barrel  eeno 

VFLPT=CHANGE  in  VflTCI'FY  due  to  change  in  propellant  1E“PERATLB  = 

xeL9NC=  PAPREL-PENO  IN  X AXIS  

XEPR0R=MISS  OISTANGE  in  X AXIS  AT  TARGET  PLANE 


XCAA IM=CUNNER  “UZZLE  AIM  IN  X AXIS 
-igiApj  FL  bend  fFEOGT  


vSIPNC^BARr Fl  RENO  ffeEGT  IN  Y AXIS 
YEPR0P  = MI.SS  OISTANCF  in  Y axis  at  target  PLANF 
YCNA  IMgCyNNFR  MUZZL'^  fl4M— 4 N-Y— AX  I S - 


DIMFNSION  ATMP (15) ,CWNC( 15 ) .pa  I N ( 15 ) .VAPRFS (15I,EFIZ,IS), 
lFTTPPC(P),FTRNG  = U),FTVLFY(AI,FTONTYUl,FTATMP(i,),FTBGHC(«il, 

2FTGRMC  (A  » iP  TCLEVIA  I ; <=  T CRFT  ( A ) t F TS  PEL  I P I . STOVEL  ( A I i AL  TO  = F ( A I , 

3P  ■'PX  I A I ,R  TR  Yl  P I ,PTBx  ( PI  ,BTBY(P)  ,OT  OX  ( P I , 0 TO  Y IP  ) , 

POLMHTJjj.NRCnj  ,NP0S(3),A|6l,Bj2)l  ,C123! 

Off-ENSICN  tFTUao)  , U-^TdB.IO) 

or  n j=i,is 

BFAC-(-5-r4-)0)  t.-A-t.M-P-U-)-^WN4U-)  .-RA-INUj  )-,VAPRF-SLjt-^F(-lW  I ,FF  13-,.J4- 


1000  fopmaT lEFin.OI 
10  CCNTINLE 
- - DC  ?C  '<  = UP^ 

P=AC(S.llOOI  FTTPPOIKI  .FTRnCEIXI  .FTVLTYIK)  .FirMYIXI  ,FTATMB(K), 
IFTPC-HC !K),FTCRwn(K),FT  = LEV(K),FTOPFT|X),FTSPFLlN)  .STOVELIK)  , 

2iL  IoefIk) 

1100  ECFMAT ( IZFft .0) 

ji).  rCAT  INUF- 


RTRYILI.RTBXIL)  .BTRYILI  .CTrx(L)  .CTCYILI 


or  30  l’l,P 
RGACIB,  12001  RTRX(L) 

1200  F'lPMAT  (BEl  JiO)- 

30  rCNTINLE 

P AT  ( S I 1300  1 ALEN,  ALPHA.RD  ,R  I , ( CUMHT  ( 1 1 , ! = I . 2 )' 

1333  FCfMaf IaFIo.O) 

PEADI'.IPOO)  STPNTY,STA’’Mp,\CT»i,n,XrNa!v,xTNFrS,Y''NFCS,V‘:iOT, 
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iv^LPetCCMi  ,cn^^,cn^^^ 

— 

P''/0  t 5 1 150n  I NnBST,(NRn';(II.I  = l,3l,^E^G 
1500  F:PMftTi5I10) 

nc  40  Kz=i  ,<i 

o = iC(5,SnCOI(lFT(KZ,K'fl,KY  = l,lOI 
_5O0Q  nOA'«  ) 

40  r-fSTiNue 

OG  50  KZ=l,l5 

-B-iie.(-5-,-90r.04U-HtHM.rX-V  I ,K vitriol- — 


50  rOFTIMP 

F AC4* (PC**4-P 1**4) 

ft''C!*P -PC/R I - - 

ALAP*(RI**3/ftL':C(ATrMP))*(((l,/3.)*(aTCMP’»*3l)«(4LF:r,(«T0PP)) 

l-((l./5,)*14TQMP**3))  + (l./'3.)) 

PCf N* 1018 .6*iLF^♦AL \H*auPF4/R 404 
Of’  50C  1=1,4 


OC — 400 — xl— I-  ,-1-5— 


rC  36  0 LL=2  ,3 
NNCQ^NPCS ILL  1 

- Cr'355  KK=l  i?  

r L^mT  sCL.RHT  i KK  ) 

r':,35c,i<T  = i ,3 

15  0 a'o?(KT)=0 

WPTT=(t,2lOC1(lFT(l,KY),KY  = l,10),(IWFT(J,KZ),KZ  = l,lO  ) , NA'P  C , CUf'Y  1 

;iOq  F4JiMA-TOtHl-T-9)0,-*F4g-P40— tABOC 1 ^1.044-,-/ 

11  OX,  <V,F4THE0  CASF ',  1044,/ 

Z1C>,'PCLN0S  PEP  nUPST  ’,12/ 

OICX;  ’40CITTVE  EFFECT  F4C  Tnp  ■,F5.2) 

ROITF  (6  ,22001 

_ _ 2200  F(;P04T(//lX,'ENn.  NO.  PURST  Nn,  OCUAC  M0.',12X, 

IFFIxEC  PIAS  (mile  1 ' , 15X,  • PAMO'OM  fkppr  (“lLj)',l>, 

2'TCTAL  FRROR  ( RFTERE  ) •/ 35X  , ' A/I  «UT)- • , 13X  , 'EL  Eva  T 1 CN'  , ICX  , 
j-l.AjlM04H'  r4-X  ,Jg-0OVA-HCN-'  ,-l-X-r-'42I  MIJTH-'  ,3Xr'EWEV44I0*i'  ,-/ 

414X,'eEAn',6X,'T0TAL',5X,'PFh4n',6X,'TCTAL',//1 
OC  300  x=1,AENG 

tF(K;GT;nGOT0210  

X“13AC  = .0. 

tPLR4C|C.  . _ _ 

0O25C 

21 0 CCFTIAUF 

4- C4  34  C C-T-O  P-P  0. 


PTEHE  = ATMP(J)*(  iq.8*MR0( 1 ) 1 
GCTC23C 

220-CCMIIxCE 

PT  EMP^ATMPI  J)tnR.a*(NR0(2)'-(CtlMMT*NPC(l)))) 

GC.TQZiO . 

230  rC4T [ALE 

rELTX  = FF( 1,  J)«(PTEMP-ATMP(JI ) 

C44-VY-«00-4-3-^  A4MP4-H-0 


XPLPNC=BCON*OELTX 
YEl  ENC  = PCCN'*nEL  OY 

250-  CrATIKL'E  -■ 

XP=XGNAIM*-X0LeNC*XTNFCS 

X t=FTGRFT  ( f ) *(  FTCP  Wn  ( I I •CWNO  ( J I *CGN3  I ♦ ( 'LTOPp  ( I )*X,M  ) 

Xl  1 = XT/ALT0EF ( I j 
C l=FTEPELt I )*FTSLFV(  ! ) 
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